1 Nephrops norvegicus FU25

1.1 Nephrops norvegicus FU25 Males

Life history parameters:

\(L_{mat}\) = 25.0 mm. \(L_{inf}^{LIT}\) = 86.0 mm. \({M/k}^{LIT}\) = 1.9.

Below, catch/landings at length per year and mean weights-at-length per year are loaded.

wal <- read.csv("25_Males_w.csv", stringsAsFactors = FALSE)
head(wal)
##   MeanLength     X1982     X1983     X1984     X1985     X1986     X1987
## 1       19.5  5.128352  5.128352  5.128352  5.128352  5.128352  5.128352
## 2       20.5  6.006350  6.006350  6.006350  6.006350  6.006350  6.006350
## 3       21.5  6.981906  6.981906  6.981906  6.981906  6.981906  6.981906
## 4       22.5  8.060560  8.060560  8.060560  8.060560  8.060560  8.060560
## 5       23.5  9.247895  9.247895  9.247895  9.247895  9.247895  9.247895
## 6       24.5 10.549534 10.549534 10.549534 10.549534 10.549534 10.549534
##       X1988     X1989     X1990     X1991     X1992     X1993     X1994
## 1  5.128352  5.128352  5.128352  5.128352  5.128352  5.128352  5.128352
## 2  6.006350  6.006350  6.006350  6.006350  6.006350  6.006350  6.006350
## 3  6.981906  6.981906  6.981906  6.981906  6.981906  6.981906  6.981906
## 4  8.060560  8.060560  8.060560  8.060560  8.060560  8.060560  8.060560
## 5  9.247895  9.247895  9.247895  9.247895  9.247895  9.247895  9.247895
## 6 10.549534 10.549534 10.549534 10.549534 10.549534 10.549534 10.549534
##       X1995     X1996     X1997     X1998     X1999     X2000     X2001
## 1  5.128352  5.128352  5.128352  5.128352  5.128352  5.128352  5.128352
## 2  6.006350  6.006350  6.006350  6.006350  6.006350  6.006350  6.006350
## 3  6.981906  6.981906  6.981906  6.981906  6.981906  6.981906  6.981906
## 4  8.060560  8.060560  8.060560  8.060560  8.060560  8.060560  8.060560
## 5  9.247895  9.247895  9.247895  9.247895  9.247895  9.247895  9.247895
## 6 10.549534 10.549534 10.549534 10.549534 10.549534 10.549534 10.549534
##       X2002     X2003     X2004     X2005     X2006     X2007     X2008
## 1  5.128352  5.128352  5.128352  5.128352  5.128352  5.128352  5.128352
## 2  6.006350  6.006350  6.006350  6.006350  6.006350  6.006350  6.006350
## 3  6.981906  6.981906  6.981906  6.981906  6.981906  6.981906  6.981906
## 4  8.060560  8.060560  8.060560  8.060560  8.060560  8.060560  8.060560
## 5  9.247895  9.247895  9.247895  9.247895  9.247895  9.247895  9.247895
## 6 10.549534 10.549534 10.549534 10.549534 10.549534 10.549534 10.549534
##       X2009     X2010     X2011     X2012     X2013     X2014     X2015
## 1  5.128352  5.128352  5.128352  5.128352  5.128352  5.128352  5.128352
## 2  6.006350  6.006350  6.006350  6.006350  6.006350  6.006350  6.006350
## 3  6.981906  6.981906  6.981906  6.981906  6.981906  6.981906  6.981906
## 4  8.060560  8.060560  8.060560  8.060560  8.060560  8.060560  8.060560
## 5  9.247895  9.247895  9.247895  9.247895  9.247895  9.247895  9.247895
## 6 10.549534 10.549534 10.549534 10.549534 10.549534 10.549534 10.549534
##       X2016     X2017     X2018     X2019
## 1  5.128352  5.128352  5.128352  5.128352
## 2  6.006350  6.006350  6.006350  6.006350
## 3  6.981906  6.981906  6.981906  6.981906
## 4  8.060560  8.060560  8.060560  8.060560
## 5  9.247895  9.247895  9.247895  9.247895
## 6 10.549534 10.549534 10.549534 10.549534
freq <- read.csv("25_Males.csv",  stringsAsFactors = FALSE)
head(freq)
##   MeanLength X1982  X1983 X1984 X1985 X1986 X1987 X1988 X1989 X1990 X1991 X1992
## 1       19.5  8148   5570  4949    NA    NA    NA    NA    NA    NA    NA    NA
## 2       20.5  1059   4154 15621     0 34637   438   301    NA    NA    NA    NA
## 3       21.5  3886    273     0     0  1064 49328   822     0    NA    NA    NA
## 4       22.5  4968  27978     0 14601   440     0 29690   715   628  1558  3748
## 5       23.5 17598  38825  9482 32234 26483  6506  3623  1624 11625 14845  5640
## 6       24.5 10263 136660 65080 50657 61017  3432   735  2552 11935 12802 75433
##   X1993 X1994 X1995 X1996 X1997 X1998 X1999 X2000 X2001 X2002 X2003 X2004 X2005
## 1    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 2    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 3    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 4   344   858   269   122   527    NA    NA    NA    NA    NA    NA    NA    NA
## 5     0    76  1125   358     0   286   204   375   467   234   310  16.7    NA
## 6  5303   651 14123  5902   366  2177   855   756    77   936   776 238.0   294
##   X2006 X2007 X2008 X2009 X2010 X2011 X2012 X2013 X2014 X2015 X2016 X2017 X2018
## 1    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 2    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 3    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 4    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 5    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 6    90 19.66    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
##   X2019
## 1    NA
## 2    NA
## 3    NA
## 4    NA
## 5    NA
## 6    NA
bin_plot(freq, 1, "mm")

On based on the previous plots of the length distribution BinWidth is fixed equal to 1.

Reference model:

lb_tableSH(freq, 1, "mm", linf = 86.0 , lmat = 25.0 , mk_ratio = 1.90 ,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1982

1.26

1.30

0.58

0.00

0.72

0.89

1983

1.10

1.22

0.56

0.00

0.67

0.89

1984

1.18

1.26

0.59

0.00

0.69

0.88

1985

1.18

1.26

0.58

0.00

0.70

0.89

1986

1.26

1.30

0.56

0.00

0.71

0.88

1987

1.30

1.30

0.59

0.01

0.74

0.90

1988

1.34

1.38

0.57

0.00

0.74

0.88

1989

1.38

1.46

0.63

0.01

0.79

0.93

1990

1.26

1.34

0.60

0.00

0.72

0.89

1991

1.14

1.26

0.55

0.00

0.67

0.87

1992

1.30

1.34

0.57

0.00

0.73

0.88

1993

1.30

1.38

0.59

0.00

0.74

0.90

1994

1.22

1.30

0.59

0.00

0.72

0.90

1995

1.22

1.34

0.58

0.00

0.72

0.90

1996

1.26

1.34

0.59

0.00

0.72

0.89

1997

1.22

1.34

0.56

0.00

0.70

0.88

1998

1.38

1.42

0.61

0.00

0.78

0.90

1999

1.22

1.34

0.63

0.01

0.73

0.91

2000

1.30

1.38

0.70

0.03

0.75

0.90

2001

1.30

1.34

0.62

0.01

0.74

0.89

2002

1.30

1.38

0.66

0.02

0.76

0.91

2003

1.42

1.50

0.71

0.03

0.83

0.95

2004

1.30

1.42

0.70

0.03

0.80

0.96

2005

1.30

1.38

0.70

0.03

0.77

0.93

2006

1.34

1.42

0.69

0.02

0.79

0.93

2007

1.34

1.46

0.68

0.02

0.79

0.94

2008

1.46

1.54

0.73

0.05

0.86

0.97

2009

1.54

1.62

0.73

0.10

0.92

1.00

2010

1.38

1.50

0.69

0.03

0.84

0.98

2011

1.54

1.62

0.67

0.02

0.84

0.91

2012

1.62

1.62

0.68

0.02

0.90

0.94

2013

1.58

1.50

0.68

0.03

0.86

0.93

2014

1.42

1.46

0.65

0.01

0.79

0.91

2015

1.38

1.42

0.62

0.00

0.78

0.91

2016

1.42

1.46

0.61

0.00

0.77

0.88

2017

1.10

1.42

0.68

0.03

0.77

1.02

2018

1.50

1.54

0.66

0.01

0.83

0.92

2019

1.42

1.50

0.60

0.00

0.78

0.90

Underestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 1, "mm", linf=86.0, lmat=25.0, mk_ratio=0.75*1.9,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1982

1.26

1.30

0.58

0.00

0.65

0.83

1983

1.10

1.22

0.56

0.00

0.61

0.83

1984

1.18

1.26

0.59

0.00

0.62

0.82

1985

1.18

1.26

0.58

0.00

0.63

0.83

1986

1.26

1.30

0.56

0.00

0.64

0.82

1987

1.30

1.30

0.59

0.00

0.67

0.85

1988

1.34

1.38

0.57

0.00

0.67

0.83

1989

1.38

1.46

0.63

0.00

0.72

0.87

1990

1.26

1.34

0.60

0.00

0.65

0.83

1991

1.14

1.26

0.55

0.00

0.60

0.81

1992

1.30

1.34

0.57

0.00

0.66

0.83

1993

1.30

1.38

0.59

0.00

0.67

0.84

1994

1.22

1.30

0.59

0.00

0.65

0.84

1995

1.22

1.34

0.58

0.00

0.65

0.84

1996

1.26

1.34

0.59

0.00

0.65

0.83

1997

1.22

1.34

0.56

0.00

0.63

0.82

1998

1.38

1.42

0.61

0.00

0.70

0.85

1999

1.22

1.34

0.63

0.01

0.65

0.85

2000

1.30

1.38

0.70

0.01

0.67

0.85

2001

1.30

1.34

0.62

0.00

0.66

0.84

2002

1.30

1.38

0.66

0.01

0.68

0.86

2003

1.42

1.50

0.71

0.01

0.75

0.90

2004

1.30

1.42

0.70

0.01

0.72

0.90

2005

1.30

1.38

0.70

0.01

0.70

0.87

2006

1.34

1.42

0.69

0.01

0.71

0.88

2007

1.34

1.46

0.68

0.01

0.71

0.88

2008

1.46

1.54

0.73

0.02

0.78

0.92

2009

1.54

1.62

0.73

0.02

0.83

0.95

2010

1.38

1.50

0.69

0.00

0.76

0.92

2011

1.54

1.62

0.67

0.00

0.76

0.87

2012

1.62

1.62

0.68

0.00

0.81

0.90

2013

1.58

1.50

0.68

0.00

0.78

0.88

2014

1.42

1.46

0.65

0.00

0.71

0.86

2015

1.38

1.42

0.62

0.00

0.70

0.86

2016

1.42

1.46

0.61

0.00

0.70

0.84

2017

1.10

1.42

0.68

0.00

0.70

0.95

2018

1.50

1.54

0.66

0.00

0.75

0.88

2019

1.42

1.50

0.60

0.00

0.71

0.85

Overestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 1, "mm", linf=86.0, lmat=25.0, mk_ratio=1.25*1.9,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1982

1.26

1.30

0.58

0.01

0.79

0.93

1983

1.10

1.22

0.56

0.01

0.74

0.94

1984

1.18

1.26

0.59

0.01

0.76

0.92

1985

1.18

1.26

0.58

0.01

0.77

0.93

1986

1.26

1.30

0.56

0.00

0.78

0.92

1987

1.30

1.30

0.59

0.01

0.82

0.94

1988

1.34

1.38

0.57

0.01

0.82

0.92

1989

1.38

1.46

0.63

0.03

0.87

0.96

1990

1.26

1.34

0.60

0.01

0.79

0.93

1991

1.14

1.26

0.55

0.00

0.73

0.92

1992

1.30

1.34

0.57

0.00

0.80

0.92

1993

1.30

1.38

0.59

0.01

0.82

0.94

1994

1.22

1.30

0.59

0.01

0.79

0.94

1995

1.22

1.34

0.58

0.01

0.79

0.94

1996

1.26

1.34

0.59

0.01

0.79

0.93

1997

1.22

1.34

0.56

0.00

0.77

0.92

1998

1.38

1.42

0.61

0.01

0.85

0.94

1999

1.22

1.34

0.63

0.02

0.80

0.95

2000

1.30

1.38

0.70

0.04

0.82

0.94

2001

1.30

1.34

0.62

0.02

0.81

0.93

2002

1.30

1.38

0.66

0.03

0.83

0.95

2003

1.42

1.50

0.71

0.07

0.91

0.98

2004

1.30

1.42

0.70

0.06

0.87

1.00

2005

1.30

1.38

0.70

0.05

0.85

0.97

2006

1.34

1.42

0.69

0.06

0.86

0.97

2007

1.34

1.46

0.68

0.04

0.87

0.98

2008

1.46

1.54

0.73

0.12

0.94

1.00

2009

1.54

1.62

0.73

0.23

1.01

1.04

2010

1.38

1.50

0.69

0.11

0.92

1.02

2011

1.54

1.62

0.67

0.06

0.92

0.94

2012

1.62

1.62

0.68

0.08

0.98

0.97

2013

1.58

1.50

0.68

0.09

0.95

0.96

2014

1.42

1.46

0.65

0.04

0.87

0.94

2015

1.38

1.42

0.62

0.02

0.85

0.94

2016

1.42

1.46

0.61

0.02

0.85

0.92

2017

1.10

1.42

0.68

0.04

0.84

1.08

2018

1.50

1.54

0.66

0.06

0.91

0.95

2019

1.42

1.50

0.60

0.01

0.86

0.93

M/k=1.5.

lb_tableSH(freq, 1, "m", linf=86.0, lmat=25.0, mk_ratio=1.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1982

1.26

1.30

0.58

0.00

0.66

0.84

1983

1.10

1.22

0.56

0.00

0.62

0.84

1984

1.18

1.26

0.59

0.00

0.63

0.83

1985

1.18

1.26

0.58

0.00

0.64

0.84

1986

1.26

1.30

0.56

0.00

0.66

0.83

1987

1.30

1.30

0.59

0.00

0.68

0.85

1988

1.34

1.38

0.57

0.00

0.68

0.84

1989

1.38

1.46

0.63

0.00

0.73

0.88

1990

1.26

1.34

0.60

0.00

0.66

0.84

1991

1.14

1.26

0.55

0.00

0.62

0.82

1992

1.30

1.34

0.57

0.00

0.67

0.84

1993

1.30

1.38

0.59

0.00

0.68

0.85

1994

1.22

1.30

0.59

0.00

0.66

0.85

1995

1.22

1.34

0.58

0.00

0.66

0.85

1996

1.26

1.34

0.59

0.00

0.66

0.84

1997

1.22

1.34

0.56

0.00

0.64

0.83

1998

1.38

1.42

0.61

0.00

0.71

0.86

1999

1.22

1.34

0.63

0.01

0.67

0.86

2000

1.30

1.38

0.70

0.02

0.69

0.86

2001

1.30

1.34

0.62

0.00

0.68

0.84

2002

1.30

1.38

0.66

0.01

0.69

0.87

2003

1.42

1.50

0.71

0.01

0.76

0.90

2004

1.30

1.42

0.70

0.01

0.73

0.91

2005

1.30

1.38

0.70

0.01

0.71

0.88

2006

1.34

1.42

0.69

0.01

0.72

0.89

2007

1.34

1.46

0.68

0.01

0.73

0.89

2008

1.46

1.54

0.73

0.02

0.79

0.93

2009

1.54

1.62

0.73

0.03

0.85

0.96

2010

1.38

1.50

0.69

0.00

0.77

0.93

2011

1.54

1.62

0.67

0.00

0.77

0.88

2012

1.62

1.62

0.68

0.00

0.82

0.91

2013

1.58

1.50

0.68

0.00

0.79

0.89

2014

1.42

1.46

0.65

0.00

0.73

0.87

2015

1.38

1.42

0.62

0.00

0.71

0.86

2016

1.42

1.46

0.61

0.00

0.71

0.84

2017

1.10

1.42

0.68

0.00

0.71

0.96

2018

1.50

1.54

0.66

0.00

0.77

0.88

2019

1.42

1.50

0.60

0.00

0.72

0.86

Underestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 1, "mm", linf=0.75*86.0, lmat=25.0, mk_ratio=1.9,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1982

1.26

1.30

0.78

0.13

0.96

0.99

1983

1.10

1.22

0.75

0.08

0.90

1.01

1984

1.18

1.26

0.79

0.11

0.92

0.99

1985

1.18

1.26

0.77

0.13

0.93

1.00

1986

1.26

1.30

0.75

0.10

0.95

0.98

1987

1.30

1.30

0.79

0.17

0.99

1.00

1988

1.34

1.38

0.76

0.13

0.99

0.98

1989

1.38

1.46

0.84

0.33

1.06

1.03

1990

1.26

1.34

0.80

0.14

0.96

0.99

1991

1.14

1.26

0.74

0.08

0.89

0.98

1992

1.30

1.34

0.77

0.13

0.97

0.98

1993

1.30

1.38

0.79

0.19

0.99

1.00

1994

1.22

1.30

0.79

0.12

0.96

1.01

1995

1.22

1.34

0.77

0.13

0.96

1.01

1996

1.26

1.34

0.79

0.13

0.96

0.99

1997

1.22

1.34

0.75

0.08

0.94

0.98

1998

1.38

1.42

0.82

0.24

1.04

1.00

1999

1.22

1.34

0.84

0.16

0.97

1.02

2000

1.30

1.38

0.93

0.16

1.00

1.00

2001

1.30

1.34

0.82

0.15

0.98

0.99

2002

1.30

1.38

0.88

0.22

1.01

1.02

2003

1.42

1.50

0.94

0.40

1.10

1.05

2004

1.30

1.42

0.93

0.32

1.06

1.07

2005

1.30

1.38

0.93

0.24

1.03

1.04

2006

1.34

1.42

0.92

0.27

1.05

1.04

2007

1.34

1.46

0.90

0.32

1.05

1.04

2008

1.46

1.54

0.98

0.47

1.15

1.07

2009

1.54

1.62

0.98

0.66

1.23

1.10

2010

1.38

1.50

0.93

0.47

1.12

1.09

2011

1.54

1.62

0.89

0.45

1.12

1.01

2012

1.62

1.62

0.91

0.57

1.19

1.04

2013

1.58

1.50

0.91

0.43

1.15

1.02

2014

1.42

1.46

0.86

0.29

1.05

1.00

2015

1.38

1.42

0.82

0.29

1.04

1.00

2016

1.42

1.46

0.81

0.22

1.03

0.98

2017

1.10

1.42

0.91

0.35

1.03

1.15

2018

1.50

1.54

0.88

0.41

1.11

1.02

2019

1.42

1.50

0.80

0.26

1.05

0.99

Overestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 1, "mm", linf=1.25*86.0, lmat=25.0, mk_ratio=1.9,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1982

1.26

1.30

0.47

0

0.58

0.80

1983

1.10

1.22

0.45

0

0.54

0.80

1984

1.18

1.26

0.47

0

0.55

0.79

1985

1.18

1.26

0.46

0

0.56

0.80

1986

1.26

1.30

0.45

0

0.57

0.79

1987

1.30

1.30

0.48

0

0.60

0.81

1988

1.34

1.38

0.46

0

0.59

0.80

1989

1.38

1.46

0.50

0

0.64

0.84

1990

1.26

1.34

0.48

0

0.58

0.80

1991

1.14

1.26

0.44

0

0.54

0.78

1992

1.30

1.34

0.46

0

0.58

0.80

1993

1.30

1.38

0.47

0

0.59

0.81

1994

1.22

1.30

0.47

0

0.57

0.81

1995

1.22

1.34

0.46

0

0.58

0.81

1996

1.26

1.34

0.48

0

0.58

0.80

1997

1.22

1.34

0.45

0

0.56

0.79

1998

1.38

1.42

0.49

0

0.62

0.82

1999

1.22

1.34

0.51

0

0.58

0.82

2000

1.30

1.38

0.56

0

0.60

0.82

2001

1.30

1.34

0.49

0

0.59

0.81

2002

1.30

1.38

0.53

0

0.60

0.83

2003

1.42

1.50

0.57

0

0.66

0.86

2004

1.30

1.42

0.56

0

0.64

0.87

2005

1.30

1.38

0.56

0

0.62

0.84

2006

1.34

1.42

0.55

0

0.63

0.85

2007

1.34

1.46

0.54

0

0.63

0.85

2008

1.46

1.54

0.59

0

0.69

0.88

2009

1.54

1.62

0.59

0

0.74

0.92

2010

1.38

1.50

0.56

0

0.67

0.89

2011

1.54

1.62

0.53

0

0.67

0.84

2012

1.62

1.62

0.54

0

0.72

0.87

2013

1.58

1.50

0.54

0

0.69

0.85

2014

1.42

1.46

0.52

0

0.63

0.82

2015

1.38

1.42

0.49

0

0.62

0.82

2016

1.42

1.46

0.49

0

0.62

0.80

2017

1.10

1.42

0.54

0

0.62

0.92

2018

1.50

1.54

0.53

0

0.67

0.84

2019

1.42

1.50

0.48

0

0.63

0.82

1.2 Nephrops norvegicus FU25 Females

Life history parameters:

\(L_{mat}\) = 28.0 mm. \(L_{inf}^{LIT}\) = 71.0 mm. \({M/k}^{LIT}\) = 1.9.

Below, catch/landings at length per year and mean weights-at-length per year are loaded.

wal <- read.csv("25_Females_w.csv", stringsAsFactors = FALSE)
freq <- read.csv("25_Females.csv",  stringsAsFactors = FALSE)
head(wal)
##   MeanLength    X1982    X1983    X1984    X1985    X1986    X1987    X1988
## 1       15.5 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 2       16.5 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 3       17.5 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
## 4       18.5 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413
## 5       19.5 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352
## 6       20.5 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350
##      X1989    X1990    X1991    X1992    X1993    X1994    X1995    X1996
## 1 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 2 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 3 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
## 4 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413
## 5 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352
## 6 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350
##      X1997    X1998    X1999    X2000    X2001    X2002    X2003    X2004
## 1 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 2 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 3 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
## 4 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413
## 5 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352
## 6 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350
##      X2005    X2006    X2007    X2008    X2009    X2010    X2011    X2012
## 1 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 2 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 3 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
## 4 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413
## 5 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352
## 6 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350
##      X2013    X2014    X2015    X2016    X2017    X2018    X2019
## 1 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 2 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 3 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
## 4 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413 4.342413
## 5 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352 5.128352
## 6 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350 6.006350
head(freq)
##   MeanLength X1982 X1983 X1984 X1985 X1986 X1987 X1988 X1989 X1990 X1991 X1992
## 1       15.5  6044    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 2       16.5     0    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 3       17.5     0    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 4       18.5     0    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 5       19.5   791     0    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 6       20.5     0 13378   782  2494   925    77   223    NA    NA    NA    NA
##   X1993 X1994 X1995 X1996 X1997 X1998 X1999 X2000 X2001 X2002 X2003 X2004 X2005
## 1    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 2    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 3    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 4    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 5    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 6    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
##   X2006 X2007 X2008 X2009 X2010 X2011 X2012 X2013 X2014 X2015 X2016 X2017 X2018
## 1    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 2    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 3    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 4    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 5    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
## 6    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
##   X2019
## 1    NA
## 2    NA
## 3    NA
## 4    NA
## 5    NA
## 6    NA
bin_plot(freq, 1, "mm")

On based on the previous plots of the length distribution BinWidth is fixed equal to 1.

Reference model:

lb_tableSH(freq, 1, "mm", linf = 71.0 , lmat = 28.0 , mk_ratio = 1.90 ,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1982

1.09

1.12

0.67

0.02

0.83

0.93

1983

0.91

1.02

0.62

0.00

0.75

0.94

1984

1.05

1.05

0.64

0.01

0.78

0.89

1985

0.98

1.05

0.65

0.01

0.77

0.91

1986

1.02

1.09

0.66

0.02

0.80

0.94

1987

0.98

1.12

0.67

0.02

0.83

0.98

1988

1.09

1.16

0.69

0.03

0.84

0.94

1989

1.23

1.27

0.69

0.03

0.91

0.94

1990

0.98

1.12

0.69

0.02

0.84

1.00

1991

1.02

1.09

0.66

0.01

0.79

0.92

1992

1.09

1.16

0.63

0.00

0.82

0.91

1993

1.20

1.20

0.65

0.01

0.86

0.91

1994

1.02

1.12

0.63

0.01

0.81

0.94

1995

1.09

1.16

0.64

0.01

0.84

0.94

1996

1.09

1.12

0.67

0.01

0.83

0.92

1997

1.09

1.16

0.65

0.01

0.82

0.92

1998

1.09

1.23

0.71

0.03

0.88

0.98

1999

1.09

1.20

0.69

0.03

0.86

0.96

2000

1.16

1.20

0.71

0.03

0.87

0.94

2001

1.05

1.16

0.67

0.01

0.84

0.95

2002

1.16

1.23

0.69

0.03

0.87

0.93

2003

1.30

1.30

0.72

0.05

0.93

0.93

2004

1.16

1.20

0.71

0.04

0.89

0.96

2005

1.12

1.20

0.71

0.04

0.88

0.97

2006

1.12

1.20

0.73

0.06

0.89

0.97

2007

1.20

1.27

0.72

0.05

0.91

0.96

2008

1.27

1.30

0.76

0.11

0.97

0.98

2009

0.91

1.38

0.77

0.22

0.98

1.22

2010

1.20

1.27

0.75

0.11

0.95

1.00

2011

1.38

1.34

0.73

0.08

0.96

0.92

2012

1.27

1.23

0.76

0.09

0.96

0.97

2013

1.02

1.16

0.73

0.07

0.87

1.01

2014

1.23

1.27

0.74

0.08

0.93

0.96

2015

1.09

1.20

0.71

0.05

0.88

0.98

2016

1.27

1.27

0.66

0.01

0.89

0.90

2017

1.23

1.30

0.76

0.09

0.94

0.97

2018

1.30

1.34

0.71

0.06

0.96

0.96

2019

1.23

1.30

0.70

0.04

0.91

0.94

Underestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 1, "mm", linf=71.0, lmat=28.0, mk_ratio=0.75*1.9,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1982

1.09

1.12

0.67

0.00

0.75

0.88

1983

0.91

1.02

0.62

0.00

0.68

0.88

1984

1.05

1.05

0.64

0.00

0.71

0.85

1985

0.98

1.05

0.65

0.00

0.69

0.86

1986

1.02

1.09

0.66

0.01

0.73

0.88

1987

0.98

1.12

0.67

0.00

0.75

0.92

1988

1.09

1.16

0.69

0.00

0.76

0.89

1989

1.23

1.27

0.69

0.00

0.82

0.90

1990

0.98

1.12

0.69

0.00

0.76

0.94

1991

1.02

1.09

0.66

0.00

0.71

0.87

1992

1.09

1.16

0.63

0.00

0.74

0.87

1993

1.20

1.20

0.65

0.00

0.78

0.87

1994

1.02

1.12

0.63

0.00

0.73

0.89

1995

1.09

1.16

0.64

0.00

0.76

0.89

1996

1.09

1.12

0.67

0.00

0.75

0.88

1997

1.09

1.16

0.65

0.00

0.74

0.87

1998

1.09

1.23

0.71

0.01

0.79

0.93

1999

1.09

1.20

0.69

0.00

0.78

0.91

2000

1.16

1.20

0.71

0.01

0.79

0.89

2001

1.05

1.16

0.67

0.00

0.75

0.90

2002

1.16

1.23

0.69

0.01

0.79

0.89

2003

1.30

1.30

0.72

0.01

0.84

0.89

2004

1.16

1.20

0.71

0.01

0.81

0.91

2005

1.12

1.20

0.71

0.01

0.80

0.92

2006

1.12

1.20

0.73

0.02

0.80

0.92

2007

1.20

1.27

0.72

0.01

0.82

0.92

2008

1.27

1.30

0.76

0.04

0.87

0.94

2009

0.91

1.38

0.77

0.06

0.88

1.14

2010

1.20

1.27

0.75

0.02

0.85

0.95

2011

1.38

1.34

0.73

0.01

0.87

0.89

2012

1.27

1.23

0.76

0.03

0.87

0.93

2013

1.02

1.16

0.73

0.02

0.79

0.96

2014

1.23

1.27

0.74

0.02

0.84

0.92

2015

1.09

1.20

0.71

0.00

0.79

0.93

2016

1.27

1.27

0.66

0.00

0.80

0.86

2017

1.23

1.30

0.76

0.06

0.84

0.92

2018

1.30

1.34

0.71

0.01

0.87

0.92

2019

1.23

1.30

0.70

0.00

0.82

0.90

Overestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 1, "mm", linf=71.0, lmat=28.0, mk_ratio=1.25*1.9,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1982

1.09

1.12

0.67

0.06

0.92

0.97

1983

0.91

1.02

0.62

0.02

0.83

0.98

1984

1.05

1.05

0.64

0.03

0.86

0.93

1985

0.98

1.05

0.65

0.03

0.84

0.95

1986

1.02

1.09

0.66

0.03

0.88

0.97

1987

0.98

1.12

0.67

0.09

0.91

1.02

1988

1.09

1.16

0.69

0.08

0.93

0.98

1989

1.23

1.27

0.69

0.10

1.00

0.97

1990

0.98

1.12

0.69

0.10

0.92

1.04

1991

1.02

1.09

0.66

0.05

0.86

0.95

1992

1.09

1.16

0.63

0.02

0.90

0.95

1993

1.20

1.20

0.65

0.04

0.95

0.94

1994

1.02

1.12

0.63

0.02

0.89

0.98

1995

1.09

1.16

0.64

0.04

0.92

0.97

1996

1.09

1.12

0.67

0.05

0.91

0.96

1997

1.09

1.16

0.65

0.03

0.90

0.96

1998

1.09

1.23

0.71

0.10

0.96

1.01

1999

1.09

1.20

0.69

0.08

0.94

0.99

2000

1.16

1.20

0.71

0.08

0.96

0.97

2001

1.05

1.16

0.67

0.06

0.92

0.99

2002

1.16

1.23

0.69

0.08

0.96

0.97

2003

1.30

1.30

0.72

0.12

1.02

0.95

2004

1.16

1.20

0.71

0.11

0.98

0.99

2005

1.12

1.20

0.71

0.12

0.97

1.00

2006

1.12

1.20

0.73

0.13

0.97

1.00

2007

1.20

1.27

0.72

0.13

1.00

0.99

2008

1.27

1.30

0.76

0.24

1.06

1.01

2009

0.91

1.38

0.77

0.46

1.07

1.27

2010

1.20

1.27

0.75

0.27

1.04

1.03

2011

1.38

1.34

0.73

0.12

1.06

0.95

2012

1.27

1.23

0.76

0.21

1.05

1.00

2013

1.02

1.16

0.73

0.15

0.95

1.05

2014

1.23

1.27

0.74

0.18

1.02

0.99

2015

1.09

1.20

0.71

0.14

0.96

1.02

2016

1.27

1.27

0.66

0.05

0.98

0.93

2017

1.23

1.30

0.76

0.16

1.03

1.00

2018

1.30

1.34

0.71

0.22

1.06

0.99

2019

1.23

1.30

0.70

0.13

1.00

0.97

M/k=1.5.

lb_tableSH(freq, 1, "m", linf=71.0, lmat=28.0, mk_ratio=1.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1982

1.09

1.12

0.67

0.00

0.77

0.89

1983

0.91

1.02

0.62

0.00

0.69

0.89

1984

1.05

1.05

0.64

0.00

0.72

0.85

1985

0.98

1.05

0.65

0.00

0.71

0.87

1986

1.02

1.09

0.66

0.01

0.74

0.89

1987

0.98

1.12

0.67

0.00

0.76

0.94

1988

1.09

1.16

0.69

0.01

0.77

0.90

1989

1.23

1.27

0.69

0.01

0.84

0.91

1990

0.98

1.12

0.69

0.01

0.77

0.95

1991

1.02

1.09

0.66

0.00

0.72

0.87

1992

1.09

1.16

0.63

0.00

0.75

0.88

1993

1.20

1.20

0.65

0.00

0.79

0.88

1994

1.02

1.12

0.63

0.00

0.74

0.90

1995

1.09

1.16

0.64

0.00

0.77

0.90

1996

1.09

1.12

0.67

0.00

0.76

0.89

1997

1.09

1.16

0.65

0.00

0.76

0.88

1998

1.09

1.23

0.71

0.01

0.80

0.94

1999

1.09

1.20

0.69

0.00

0.79

0.92

2000

1.16

1.20

0.71

0.01

0.80

0.90

2001

1.05

1.16

0.67

0.00

0.77

0.91

2002

1.16

1.23

0.69

0.01

0.80

0.90

2003

1.30

1.30

0.72

0.01

0.86

0.90

2004

1.16

1.20

0.71

0.01

0.82

0.92

2005

1.12

1.20

0.71

0.01

0.81

0.93

2006

1.12

1.20

0.73

0.02

0.81

0.93

2007

1.20

1.27

0.72

0.01

0.84

0.92

2008

1.27

1.30

0.76

0.05

0.89

0.95

2009

0.91

1.38

0.77

0.08

0.90

1.16

2010

1.20

1.27

0.75

0.04

0.87

0.96

2011

1.38

1.34

0.73

0.01

0.88

0.90

2012

1.27

1.23

0.76

0.03

0.88

0.94

2013

1.02

1.16

0.73

0.02

0.80

0.97

2014

1.23

1.27

0.74

0.03

0.85

0.92

2015

1.09

1.20

0.71

0.01

0.81

0.94

2016

1.27

1.27

0.66

0.00

0.82

0.87

2017

1.23

1.30

0.76

0.06

0.86

0.93

2018

1.30

1.34

0.71

0.01

0.88

0.93

2019

1.23

1.30

0.70

0.01

0.84

0.91

Underestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 1, "mm", linf=0.75*71.0, lmat=28.0, mk_ratio=1.9,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1982

1.09

1.12

0.89

0.40

1.11

1.03

1983

0.91

1.02

0.83

0.23

1.00

1.05

1984

1.05

1.05

0.85

0.19

1.05

0.99

1985

0.98

1.05

0.87

0.21

1.02

1.02

1986

1.02

1.09

0.88

0.31

1.07

1.04

1987

0.98

1.12

0.90

0.42

1.10

1.09

1988

1.09

1.16

0.92

0.43

1.12

1.04

1989

1.23

1.27

0.93

0.68

1.21

1.03

1990

0.98

1.12

0.91

0.48

1.12

1.11

1991

1.02

1.09

0.88

0.28

1.05

1.02

1992

1.09

1.16

0.84

0.36

1.09

1.01

1993

1.20

1.20

0.86

0.49

1.15

1.00

1994

1.02

1.12

0.84

0.37

1.08

1.05

1995

1.09

1.16

0.86

0.44

1.12

1.04

1996

1.09

1.12

0.89

0.35

1.10

1.02

1997

1.09

1.16

0.87

0.38

1.10

1.02

1998

1.09

1.23

0.94

0.59

1.17

1.08

1999

1.09

1.20

0.92

0.53

1.15

1.06

2000

1.16

1.20

0.94

0.49

1.16

1.03

2001

1.05

1.16

0.89

0.44

1.11

1.05

2002

1.16

1.23

0.92

0.53

1.16

1.03

2003

1.30

1.30

0.96

0.76

1.24

1.01

2004

1.16

1.20

0.94

0.60

1.19

1.06

2005

1.12

1.20

0.95

0.57

1.18

1.07

2006

1.12

1.20

0.98

0.58

1.18

1.07

2007

1.20

1.27

0.96

0.69

1.22

1.05

2008

1.27

1.30

1.02

0.80

1.29

1.07

2009

0.91

1.38

1.02

0.86

1.31

1.36

2010

1.20

1.27

1.00

0.67

1.26

1.09

2011

1.38

1.34

0.97

0.83

1.28

1.01

2012

1.27

1.23

1.02

0.67

1.28

1.06

2013

1.02

1.16

0.98

0.52

1.16

1.12

2014

1.23

1.27

0.99

0.74

1.24

1.05

2015

1.09

1.20

0.95

0.57

1.17

1.09

2016

1.27

1.27

0.88

0.63

1.19

0.99

2017

1.23

1.30

1.02

0.77

1.25

1.06

2018

1.30

1.34

0.95

0.87

1.28

1.05

2019

1.23

1.30

0.93

0.76

1.22

1.03

Overestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 1, "mm", linf=1.25*71.0, lmat=28.0, mk_ratio=1.9,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1982

1.09

1.12

0.54

0

0.67

0.85

1983

0.91

1.02

0.50

0

0.60

0.85

1984

1.05

1.05

0.51

0

0.63

0.81

1985

0.98

1.05

0.52

0

0.61

0.83

1986

1.02

1.09

0.53

0

0.64

0.85

1987

0.98

1.12

0.54

0

0.66

0.89

1988

1.09

1.16

0.55

0

0.67

0.86

1989

1.23

1.27

0.56

0

0.73

0.86

1990

0.98

1.12

0.55

0

0.67

0.91

1991

1.02

1.09

0.53

0

0.63

0.83

1992

1.09

1.16

0.50

0

0.65

0.83

1993

1.20

1.20

0.52

0

0.69

0.83

1994

1.02

1.12

0.50

0

0.65

0.86

1995

1.09

1.16

0.51

0

0.67

0.86

1996

1.09

1.12

0.53

0

0.66

0.84

1997

1.09

1.16

0.52

0

0.66

0.84

1998

1.09

1.23

0.57

0

0.70

0.89

1999

1.09

1.20

0.55

0

0.69

0.88

2000

1.16

1.20

0.56

0

0.70

0.86

2001

1.05

1.16

0.54

0

0.67

0.87

2002

1.16

1.23

0.55

0

0.70

0.86

2003

1.30

1.30

0.57

0

0.75

0.86

2004

1.16

1.20

0.57

0

0.72

0.88

2005

1.12

1.20

0.57

0

0.71

0.88

2006

1.12

1.20

0.59

0

0.71

0.89

2007

1.20

1.27

0.57

0

0.73

0.88

2008

1.27

1.30

0.61

0

0.77

0.90

2009

0.91

1.38

0.61

0

0.78

1.10

2010

1.20

1.27

0.60

0

0.76

0.91

2011

1.38

1.34

0.58

0

0.77

0.86

2012

1.27

1.23

0.61

0

0.77

0.90

2013

1.02

1.16

0.59

0

0.70

0.92

2014

1.23

1.27

0.59

0

0.74

0.88

2015

1.09

1.20

0.57

0

0.70

0.90

2016

1.27

1.27

0.53

0

0.71

0.83

2017

1.23

1.30

0.61

0

0.75

0.89

2018

1.30

1.34

0.57

0

0.77

0.88

2019

1.23

1.30

0.56

0

0.73

0.87

2 Nephrops norvegicus FU2627

2.1 Nephrops norvegicus FU2627 Males

Life history parameters:

\(L_{mat}\) = 28.0 mm. \(L_{inf}^{LIT}\) = 80.0 mm. \({M/k}^{LIT}\) = 2.

Below, catch/landings at length per year and mean weights-at-length per year are loaded.

wal <- read.csv("FU2627_PesosMales_1988_2019_medio.csv", stringsAsFactors = FALSE)
freq <- read.csv("FU2627_tallasMales_1988_2019.csv", stringsAsFactors = FALSE)
head(freq)
##   talla X1988  X1989  X1990 X1991 X1992 X1993 X1994 X1995 X1996 X1997 X1998
## 1  12.5     0      0   7594     0     0     0     0     0     0     0     0
## 2  13.5     0  66181   8178     0     0     0     0     0     0     0     0
## 3  14.5     0  66181  10497     0     0     0     0     0     0     0     0
## 4  15.5     0 339166  84532     0     0     0     0     0     0     0     0
## 5  16.5     0  75996 135341  5393     0     0     0     0     0     0     0
## 6  17.5     0  92508 175731  8053     0     0  5015     0     0     0     0
##   X1999 X2000 X2001 X2002 X2003 X2004 X2005 X2006 X2007 X2008 X2009 X2010 X2011
## 1     0     0     0     0     0     0     0     0     0     0     0     0     0
## 2     0     0     0     0     0     0     0     0     0     0     0     0     0
## 3     0     0     0     0     0     0     0     0     0     0     0     0     0
## 4     0     0     0     0     0     0     0     0     0     0     0     0     0
## 5     0     0  4030     0     0     0     0     0     0     0     0     0     0
## 6     0     0   232   142     0     0     0     0     0     0     0     0     0
##   X2012 X2013 X2014 X2015 X2016 X2017 X2018 X2019
## 1     0     0     0     0     0     0     0     0
## 2     0     0     0     0     0     0     0     0
## 3     0     0     0     0     0     0     0     0
## 4     0     0     0     0     0     0     0     0
## 5     0     0     0     0     0     0     0     0
## 6     0     0     0     0     0     0     0     0
head(wal)
##   talla    X1988    X1989    X1990    X1991    X1992    X1993    X1994    X1995
## 1  12.5 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068
## 2  13.5 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439
## 3  14.5 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904
## 4  15.5 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 5  16.5 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 6  17.5 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
##      X1996    X1997    X1998    X1999    X2000    X2001    X2002    X2003
## 1 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068
## 2 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439
## 3 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904
## 4 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 5 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 6 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
##      X2004    X2005    X2006    X2007    X2008    X2009    X2010    X2011
## 1 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068
## 2 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439
## 3 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904
## 4 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 5 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 6 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
##      X2012    X2013    X2014    X2015    X2016    X2017    X2018    X2019
## 1 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068
## 2 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439
## 3 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904
## 4 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 5 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 6 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
bin_plot(freq, 1, "mm")

bin_plot(freq, 2, "mm")

On based on the previous plots of the length distribution BinWidth is fixed equal to 2.

Reference model:

lb_tableSH(freq, 2, "mm", linf=80, lmat=28, mk_ratio=2,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

1.11

1.11

0.62

0.01

0.75

0.89

1989

0.82

0.89

0.62

0.01

0.65

0.91

1990

0.75

0.82

0.53

0.00

0.56

0.82

1991

0.89

0.96

0.59

0.01

0.67

0.90

1992

0.96

1.11

0.65

0.02

0.77

0.98

1993

0.96

1.04

0.63

0.01

0.72

0.92

1994

1.04

1.11

0.69

0.03

0.77

0.94

1995

1.04

1.04

0.58

0.00

0.72

0.88

1996

0.89

0.96

0.61

0.01

0.68

0.90

1997

1.04

1.11

0.65

0.02

0.78

0.96

1998

1.11

1.18

0.67

0.03

0.82

0.96

1999

1.11

1.18

0.67

0.03

0.81

0.95

2000

0.82

1.04

0.73

0.04

0.74

1.03

2001

0.82

0.89

0.64

0.02

0.64

0.90

2002

1.11

1.18

0.75

0.06

0.83

0.98

2003

1.11

1.18

0.79

0.08

0.86

1.01

2004

1.25

1.32

0.79

0.09

0.90

0.99

2005

1.11

1.25

0.87

0.17

0.92

1.08

2006

1.18

1.32

0.85

0.11

0.89

1.01

2007

1.18

1.32

0.81

0.12

0.91

1.03

2008

1.18

1.32

0.84

0.19

0.94

1.07

2009

1.18

1.32

0.87

0.29

0.97

1.10

2010

1.54

1.61

0.91

0.48

1.14

1.09

2011

1.39

1.46

0.91

0.31

1.05

1.07

2012

1.39

1.54

0.91

0.31

1.06

1.08

2013

1.11

1.32

0.89

0.37

1.00

1.18

2014

1.61

1.46

0.88

0.30

1.12

1.03

2015

1.25

1.39

0.88

0.22

0.95

1.04

2016

1.25

1.39

0.86

0.28

0.99

1.09

2017

1.39

1.46

0.90

0.36

1.07

1.09

2018

1.25

1.39

0.86

0.35

1.04

1.13

2019

1.25

1.46

0.86

0.27

1.00

1.09

Underestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 2, "mm", linf=80, lmat=28, mk_ratio=0.75*2,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

1.11

1.11

0.62

0.00

0.68

0.84

1989

0.82

0.89

0.62

0.00

0.59

0.84

1990

0.75

0.82

0.53

0.00

0.50

0.75

1991

0.89

0.96

0.59

0.00

0.61

0.83

1992

0.96

1.11

0.65

0.00

0.69

0.92

1993

0.96

1.04

0.63

0.00

0.65

0.86

1994

1.04

1.11

0.69

0.01

0.69

0.89

1995

1.04

1.04

0.58

0.00

0.65

0.83

1996

0.89

0.96

0.61

0.00

0.61

0.84

1997

1.04

1.11

0.65

0.00

0.70

0.90

1998

1.11

1.18

0.67

0.01

0.73

0.91

1999

1.11

1.18

0.67

0.00

0.73

0.90

2000

0.82

1.04

0.73

0.02

0.67

0.96

2001

0.82

0.89

0.64

0.01

0.58

0.83

2002

1.11

1.18

0.75

0.03

0.75

0.92

2003

1.11

1.18

0.79

0.03

0.77

0.95

2004

1.25

1.32

0.79

0.04

0.81

0.94

2005

1.11

1.25

0.87

0.09

0.83

1.02

2006

1.18

1.32

0.85

0.06

0.80

0.96

2007

1.18

1.32

0.81

0.05

0.82

0.97

2008

1.18

1.32

0.84

0.09

0.85

1.01

2009

1.18

1.32

0.87

0.15

0.88

1.04

2010

1.54

1.61

0.91

0.25

1.03

1.05

2011

1.39

1.46

0.91

0.17

0.95

1.03

2012

1.39

1.54

0.91

0.18

0.95

1.03

2013

1.11

1.32

0.89

0.31

0.90

1.11

2014

1.61

1.46

0.88

0.20

1.01

1.00

2015

1.25

1.39

0.88

0.12

0.86

0.99

2016

1.25

1.39

0.86

0.16

0.90

1.03

2017

1.39

1.46

0.90

0.21

0.96

1.04

2018

1.25

1.39

0.86

0.25

0.93

1.08

2019

1.25

1.46

0.86

0.15

0.90

1.03

Overestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 2, "mm", linf=80, lmat=28, mk_ratio=1.25*2,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

1.11

1.11

0.62

0.02

0.83

0.92

1989

0.82

0.89

0.62

0.03

0.72

0.97

1990

0.75

0.82

0.53

0.01

0.62

0.87

1991

0.89

0.96

0.59

0.01

0.74

0.94

1992

0.96

1.11

0.65

0.06

0.84

1.03

1993

0.96

1.04

0.63

0.03

0.80

0.97

1994

1.04

1.11

0.69

0.06

0.85

0.99

1995

1.04

1.04

0.58

0.01

0.79

0.92

1996

0.89

0.96

0.61

0.02

0.75

0.95

1997

1.04

1.11

0.65

0.05

0.86

1.00

1998

1.11

1.18

0.67

0.10

0.90

1.00

1999

1.11

1.18

0.67

0.08

0.89

0.99

2000

0.82

1.04

0.73

0.07

0.82

1.10

2001

0.82

0.89

0.64

0.03

0.71

0.95

2002

1.11

1.18

0.75

0.10

0.91

1.02

2003

1.11

1.18

0.79

0.15

0.94

1.05

2004

1.25

1.32

0.79

0.15

0.99

1.02

2005

1.11

1.25

0.87

0.29

1.01

1.13

2006

1.18

1.32

0.85

0.18

0.98

1.05

2007

1.18

1.32

0.81

0.21

1.00

1.07

2008

1.18

1.32

0.84

0.30

1.04

1.11

2009

1.18

1.32

0.87

0.39

1.07

1.15

2010

1.54

1.61

0.91

0.63

1.26

1.12

2011

1.39

1.46

0.91

0.44

1.16

1.10

2012

1.39

1.54

0.91

0.48

1.17

1.11

2013

1.11

1.32

0.89

0.45

1.10

1.23

2014

1.61

1.46

0.88

0.48

1.23

1.06

2015

1.25

1.39

0.88

0.28

1.05

1.08

2016

1.25

1.39

0.86

0.40

1.09

1.12

2017

1.39

1.46

0.90

0.50

1.18

1.12

2018

1.25

1.39

0.86

0.44

1.14

1.17

2019

1.25

1.46

0.86

0.42

1.09

1.12

M/k=1.5.

lb_tableSH(freq, 2, "mm", linf=80, lmat=28, mk_ratio=1.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

1.11

1.11

0.62

0.00

0.68

0.84

1989

0.82

0.89

0.62

0.00

0.59

0.84

1990

0.75

0.82

0.53

0.00

0.50

0.75

1991

0.89

0.96

0.59

0.00

0.61

0.83

1992

0.96

1.11

0.65

0.00

0.69

0.92

1993

0.96

1.04

0.63

0.00

0.65

0.86

1994

1.04

1.11

0.69

0.01

0.69

0.89

1995

1.04

1.04

0.58

0.00

0.65

0.83

1996

0.89

0.96

0.61

0.00

0.61

0.84

1997

1.04

1.11

0.65

0.00

0.70

0.90

1998

1.11

1.18

0.67

0.01

0.73

0.91

1999

1.11

1.18

0.67

0.00

0.73

0.90

2000

0.82

1.04

0.73

0.02

0.67

0.96

2001

0.82

0.89

0.64

0.01

0.58

0.83

2002

1.11

1.18

0.75

0.03

0.75

0.92

2003

1.11

1.18

0.79

0.03

0.77

0.95

2004

1.25

1.32

0.79

0.04

0.81

0.94

2005

1.11

1.25

0.87

0.09

0.83

1.02

2006

1.18

1.32

0.85

0.06

0.80

0.96

2007

1.18

1.32

0.81

0.05

0.82

0.97

2008

1.18

1.32

0.84

0.09

0.85

1.01

2009

1.18

1.32

0.87

0.15

0.88

1.04

2010

1.54

1.61

0.91

0.25

1.03

1.05

2011

1.39

1.46

0.91

0.17

0.95

1.03

2012

1.39

1.54

0.91

0.18

0.95

1.03

2013

1.11

1.32

0.89

0.31

0.90

1.11

2014

1.61

1.46

0.88

0.20

1.01

1.00

2015

1.25

1.39

0.88

0.12

0.86

0.99

2016

1.25

1.39

0.86

0.16

0.90

1.03

2017

1.39

1.46

0.90

0.21

0.96

1.04

2018

1.25

1.39

0.86

0.25

0.93

1.08

2019

1.25

1.46

0.86

0.15

0.90

1.03

Underestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 2, "mm", linf=0.75*80, lmat=28, mk_ratio=2,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

1.11

1.11

0.82

0.14

1.00

0.98

1989

0.82

0.89

0.83

0.11

0.87

1.03

1990

0.75

0.82

0.71

0.03

0.75

0.93

1991

0.89

0.96

0.78

0.09

0.90

1.01

1992

0.96

1.11

0.87

0.30

1.02

1.10

1993

0.96

1.04

0.84

0.18

0.97

1.03

1994

1.04

1.11

0.91

0.24

1.03

1.05

1995

1.04

1.04

0.78

0.12

0.96

0.98

1996

0.89

0.96

0.81

0.13

0.90

1.02

1997

1.04

1.11

0.87

0.31

1.04

1.06

1998

1.11

1.18

0.90

0.36

1.09

1.06

1999

1.11

1.18

0.89

0.32

1.08

1.05

2000

0.82

1.04

0.98

0.23

0.99

1.17

2001

0.82

0.89

0.86

0.10

0.86

1.02

2002

1.11

1.18

1.00

0.36

1.11

1.08

2003

1.11

1.18

1.05

0.39

1.14

1.12

2004

1.25

1.32

1.06

0.55

1.21

1.09

2005

1.11

1.25

1.16

0.56

1.22

1.20

2006

1.18

1.32

1.13

0.52

1.19

1.12

2007

1.18

1.32

1.08

0.62

1.21

1.13

2008

1.18

1.32

1.12

0.63

1.26

1.18

2009

1.18

1.32

1.16

0.67

1.30

1.22

2010

1.54

1.61

1.21

0.86

1.53

1.18

2011

1.39

1.46

1.21

0.81

1.40

1.17

2012

1.39

1.54

1.21

0.82

1.41

1.18

2013

1.11

1.32

1.19

0.61

1.34

1.31

2014

1.61

1.46

1.17

0.77

1.49

1.12

2015

1.25

1.39

1.17

0.69

1.27

1.14

2016

1.25

1.39

1.15

0.74

1.33

1.19

2017

1.39

1.46

1.20

0.84

1.43

1.19

2018

1.25

1.39

1.15

0.72

1.38

1.25

2019

1.25

1.46

1.14

0.76

1.33

1.19

Overestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 2, "mm", linf=1.25*80, lmat=28, mk_ratio=2,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

1.11

1.11

0.49

0.00

0.60

0.81

1989

0.82

0.89

0.50

0.00

0.52

0.82

1990

0.75

0.82

0.42

0.00

0.45

0.73

1991

0.89

0.96

0.47

0.00

0.54

0.81

1992

0.96

1.11

0.52

0.00

0.61

0.89

1993

0.96

1.04

0.51

0.00

0.58

0.84

1994

1.04

1.11

0.55

0.00

0.62

0.86

1995

1.04

1.04

0.47

0.00

0.58

0.80

1996

0.89

0.96

0.49

0.00

0.54

0.81

1997

1.04

1.11

0.52

0.00

0.62

0.87

1998

1.11

1.18

0.54

0.00

0.65

0.87

1999

1.11

1.18

0.54

0.00

0.65

0.86

2000

0.82

1.04

0.59

0.01

0.59

0.93

2001

0.82

0.89

0.52

0.00

0.52

0.81

2002

1.11

1.18

0.60

0.01

0.66

0.89

2003

1.11

1.18

0.63

0.02

0.69

0.92

2004

1.25

1.32

0.64

0.01

0.72

0.90

2005

1.11

1.25

0.70

0.04

0.73

0.98

2006

1.18

1.32

0.68

0.03

0.71

0.92

2007

1.18

1.32

0.65

0.02

0.73

0.94

2008

1.18

1.32

0.67

0.03

0.75

0.98

2009

1.18

1.32

0.70

0.04

0.78

1.01

2010

1.54

1.61

0.73

0.10

0.92

1.01

2011

1.39

1.46

0.73

0.07

0.84

0.99

2012

1.39

1.54

0.73

0.10

0.85

0.99

2013

1.11

1.32

0.71

0.08

0.80

1.07

2014

1.61

1.46

0.70

0.04

0.89

0.96

2015

1.25

1.39

0.70

0.04

0.76

0.95

2016

1.25

1.39

0.69

0.05

0.80

0.99

2017

1.39

1.46

0.72

0.09

0.86

1.00

2018

1.25

1.39

0.69

0.08

0.83

1.04

2019

1.25

1.46

0.69

0.04

0.80

1.00

2.2 Nephrops norvegicus FU2627 Females

Life history parameters:

\(L_{mat}\) = 26.0 mm. \(L_{inf}^{LIT}\) = 65.0 mm. \({M/k}^{LIT}\) = 2.5

Below, catch/landings at length per year and mean weights-at-length per year are loaded.

wal <- read.csv("FU2627_PesosFemales_1988_2019_medio.csv", stringsAsFactors = FALSE)
freq <- read.csv("FU2627_tallasFemales_1988_2019.csv", stringsAsFactors = FALSE)
head(freq)
##   talla X1988  X1989  X1990 X1991 X1992 X1993 X1994 X1995 X1996 X1997 X1998
## 1  12.5     0      0      0     0     0     0     0     0     0     0     0
## 2  13.5     0   4239  15172     0     0     0     0     0     0     0     0
## 3  14.5     0   2474  25480     0     0     0     0     0     0     0     0
## 4  15.5     0 108363  61013  1795     0     0     0     0     0     0     0
## 5  16.5     0 113495 254545  9073     0     0     0     0     0     0     0
## 6  17.5     0  34688 529385 11820     0     0  2730     0     0     0     0
##   X1999 X2000 X2001 X2002 X2003 X2004 X2005 X2006 X2007 X2008 X2009 X2010 X2011
## 1     0     0     0     0     0     0     0     0     0     0     0     0     0
## 2     0     0     0     0     0     0     0     0     0     0     0     0     0
## 3     0     0     0     0     0     0     0     0     0     0     0     0     0
## 4     0     0     0     0     0     0     0     0     0     0     0     0     0
## 5     0  3063 13552     0     0     0     0     0     0     0     0     0     0
## 6     0  3063 11684     0     0     0     0     0     0     0     0     0     0
##   X2012 X2013 X2014 X2015 X2016 X2017 X2018 X2019
## 1     0     0     0     0     0     0     0     0
## 2     0     0     0     0     0     0     0     0
## 3     0     0     0     0     0     0     0     0
## 4     0     0     0     0     0     0     0     0
## 5     0     0     0     0     0     0     0     0
## 6     0     0     0     0     0     0     0     0
head(wal)
##   talla    X1988    X1989    X1990    X1991    X1992    X1993    X1994    X1995
## 1  12.5 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068
## 2  13.5 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439
## 3  14.5 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904
## 4  15.5 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 5  16.5 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 6  17.5 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
##      X1996    X1997    X1998    X1999    X2000    X2001    X2002    X2003
## 1 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068
## 2 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439
## 3 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904
## 4 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 5 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 6 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
##      X2004    X2005    X2006    X2007    X2008    X2009    X2010    X2011
## 1 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068
## 2 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439
## 3 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904
## 4 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 5 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 6 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
##      X2012    X2013    X2014    X2015    X2016    X2017    X2018    X2019
## 1 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068 1.258068
## 2 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439 1.604439
## 3 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904 2.010904
## 4 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656 2.482656
## 5 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947 3.024947
## 6 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080 3.643080
bin_plot(freq, 1, "cm")

bin_plot(freq, 2, "cm")

On based on the previous plots of the length distribution BinWidth is fixed equal to 2.

Reference model:

lb_tableSH(freq, 2, "mm", linf=65, lmat=26, mk_ratio=2.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

1.12

1.12

0.68

0.14

0.96

0.97

1989

0.81

0.88

0.68

0.08

0.82

1.03

1990

0.81

0.81

0.60

0.02

0.73

0.91

1991

0.96

1.04

0.69

0.09

0.87

0.98

1992

1.04

1.12

0.72

0.16

0.96

1.02

1993

1.04

1.12

0.71

0.18

0.95

1.01

1994

1.12

1.19

0.74

0.25

1.00

1.02

1995

1.04

1.12

0.68

0.11

0.93

0.98

1996

0.96

1.04

0.71

0.14

0.90

1.01

1997

1.19

1.19

0.77

0.29

1.04

1.00

1998

1.19

1.27

0.80

0.43

1.09

1.05

1999

1.19

1.27

0.79

0.47

1.10

1.06

2000

0.81

1.04

0.73

0.21

0.91

1.14

2001

0.81

0.88

0.70

0.09

0.81

1.01

2002

1.04

1.19

0.81

0.39

1.04

1.11

2003

1.19

1.19

0.76

0.29

1.05

1.01

2004

1.27

1.35

0.81

0.58

1.15

1.06

2005

1.12

1.27

0.85

0.51

1.11

1.12

2006

1.27

1.35

0.83

0.52

1.12

1.04

2007

1.27

1.35

0.83

0.45

1.11

1.02

2008

1.27

1.27

0.78

0.45

1.11

1.02

2009

1.12

1.35

0.86

0.52

1.10

1.12

2010

1.42

1.50

0.90

0.80

1.26

1.08

2011

1.35

1.50

0.90

0.78

1.25

1.11

2012

1.50

1.50

0.89

0.76

1.26

1.03

2013

1.35

1.42

0.94

0.74

1.28

1.14

2014

1.19

1.35

0.89

0.60

1.18

1.14

2015

1.27

1.35

0.86

0.62

1.15

1.06

2016

1.19

1.27

0.85

0.53

1.13

1.09

2017

1.42

1.50

1.01

0.83

1.32

1.12

2018

1.19

1.35

0.86

0.51

1.13

1.09

2019

1.19

1.42

0.88

0.71

1.20

1.16

Underestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 2, "mm", linf=65, lmat=26, mk_ratio=0.75*2.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

1.12

1.12

0.68

0.02

0.85

0.93

1989

0.81

0.88

0.68

0.02

0.73

0.96

1990

0.81

0.81

0.60

0.01

0.65

0.86

1991

0.96

1.04

0.69

0.02

0.77

0.93

1992

1.04

1.12

0.72

0.05

0.85

0.97

1993

1.04

1.12

0.71

0.04

0.85

0.97

1994

1.12

1.19

0.74

0.05

0.89

0.97

1995

1.04

1.12

0.68

0.02

0.82

0.94

1996

0.96

1.04

0.71

0.03

0.80

0.95

1997

1.19

1.19

0.77

0.10

0.92

0.96

1998

1.19

1.27

0.80

0.16

0.96

1.01

1999

1.19

1.27

0.79

0.18

0.97

1.02

2000

0.81

1.04

0.73

0.05

0.81

1.07

2001

0.81

0.88

0.70

0.03

0.71

0.94

2002

1.04

1.19

0.81

0.14

0.92

1.05

2003

1.19

1.19

0.76

0.08

0.93

0.97

2004

1.27

1.35

0.81

0.22

1.02

1.02

2005

1.12

1.27

0.85

0.25

0.98

1.07

2006

1.27

1.35

0.83

0.16

1.00

1.00

2007

1.27

1.35

0.83

0.16

0.98

0.99

2008

1.27

1.27

0.78

0.16

0.98

0.99

2009

1.12

1.35

0.86

0.22

0.98

1.07

2010

1.42

1.50

0.90

0.47

1.12

1.04

2011

1.35

1.50

0.90

0.46

1.11

1.07

2012

1.50

1.50

0.89

0.31

1.12

1.01

2013

1.35

1.42

0.94

0.39

1.14

1.10

2014

1.19

1.35

0.89

0.37

1.05

1.10

2015

1.27

1.35

0.86

0.21

1.02

1.02

2016

1.19

1.27

0.85

0.30

1.00

1.05

2017

1.42

1.50

1.01

0.53

1.17

1.09

2018

1.19

1.35

0.86

0.29

1.00

1.05

2019

1.19

1.42

0.88

0.43

1.06

1.11

Overestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 2, "mm", linf=65, lmat=26, mk_ratio=1.25*2.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

1.12

1.12

0.68

0.23

1.07

1.00

1989

0.81

0.88

0.68

0.12

0.91

1.08

1990

0.81

0.81

0.60

0.04

0.81

0.96

1991

0.96

1.04

0.69

0.13

0.97

1.01

1992

1.04

1.12

0.72

0.24

1.06

1.05

1993

1.04

1.12

0.71

0.27

1.06

1.05

1994

1.12

1.19

0.74

0.35

1.12

1.05

1995

1.04

1.12

0.68

0.18

1.03

1.02

1996

0.96

1.04

0.71

0.21

1.00

1.04

1997

1.19

1.19

0.77

0.38

1.15

1.03

1998

1.19

1.27

0.80

0.56

1.21

1.08

1999

1.19

1.27

0.79

0.61

1.22

1.09

2000

0.81

1.04

0.73

0.29

1.02

1.20

2001

0.81

0.88

0.70

0.13

0.90

1.06

2002

1.04

1.19

0.81

0.48

1.16

1.14

2003

1.19

1.19

0.76

0.44

1.16

1.04

2004

1.27

1.35

0.81

0.70

1.28

1.09

2005

1.12

1.27

0.85

0.60

1.23

1.16

2006

1.27

1.35

0.83

0.65

1.25

1.07

2007

1.27

1.35

0.83

0.60

1.23

1.05

2008

1.27

1.27

0.78

0.58

1.23

1.05

2009

1.12

1.35

0.86

0.63

1.23

1.15

2010

1.42

1.50

0.90

0.88

1.41

1.10

2011

1.35

1.50

0.90

0.84

1.39

1.13

2012

1.50

1.50

0.89

0.82

1.41

1.05

2013

1.35

1.42

0.94

0.79

1.43

1.16

2014

1.19

1.35

0.89

0.66

1.32

1.18

2015

1.27

1.35

0.86

0.74

1.28

1.09

2016

1.19

1.27

0.85

0.60

1.26

1.12

2017

1.42

1.50

1.01

0.90

1.47

1.14

2018

1.19

1.35

0.86

0.60

1.26

1.12

2019

1.19

1.42

0.88

0.81

1.33

1.19

M/k=1.5.

lb_tableSH(freq, 2, "mm", linf=65, lmat=26, mk_ratio=1.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

1.12

1.12

0.68

0.01

0.78

0.89

1989

0.81

0.88

0.68

0.01

0.67

0.91

1990

0.81

0.81

0.60

0.00

0.60

0.81

1991

0.96

1.04

0.69

0.01

0.71

0.88

1992

1.04

1.12

0.72

0.01

0.78

0.93

1993

1.04

1.12

0.71

0.01

0.78

0.93

1994

1.12

1.19

0.74

0.02

0.82

0.94

1995

1.04

1.12

0.68

0.00

0.76

0.90

1996

0.96

1.04

0.71

0.01

0.74

0.91

1997

1.19

1.19

0.77

0.03

0.85

0.93

1998

1.19

1.27

0.80

0.06

0.89

0.97

1999

1.19

1.27

0.79

0.06

0.90

0.98

2000

0.81

1.04

0.73

0.02

0.75

1.01

2001

0.81

0.88

0.70

0.01

0.66

0.89

2002

1.04

1.19

0.81

0.06

0.85

1.01

2003

1.19

1.19

0.76

0.02

0.86

0.94

2004

1.27

1.35

0.81

0.07

0.94

0.99

2005

1.12

1.27

0.85

0.11

0.91

1.03

2006

1.27

1.35

0.83

0.08

0.92

0.97

2007

1.27

1.35

0.83

0.06

0.90

0.96

2008

1.27

1.27

0.78

0.04

0.90

0.96

2009

1.12

1.35

0.86

0.11

0.90

1.03

2010

1.42

1.50

0.90

0.21

1.03

1.02

2011

1.35

1.50

0.90

0.25

1.02

1.04

2012

1.50

1.50

0.89

0.20

1.03

0.99

2013

1.35

1.42

0.94

0.32

1.05

1.07

2014

1.19

1.35

0.89

0.22

0.97

1.06

2015

1.27

1.35

0.86

0.09

0.94

0.99

2016

1.19

1.27

0.85

0.16

0.93

1.01

2017

1.42

1.50

1.01

0.32

1.08

1.06

2018

1.19

1.35

0.86

0.12

0.93

1.02

2019

1.19

1.42

0.88

0.19

0.98

1.07

Underestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 2, "mm", linf=0.75*65, lmat=26, mk_ratio=2.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

1.12

1.12

0.91

0.70

1.28

1.05

1989

0.81

0.88

0.90

0.37

1.10

1.14

1990

0.81

0.81

0.79

0.14

0.97

1.01

1991

0.96

1.04

0.92

0.43

1.16

1.07

1992

1.04

1.12

0.96

0.72

1.27

1.11

1993

1.04

1.12

0.95

0.68

1.27

1.10

1994

1.12

1.19

0.98

0.79

1.34

1.10

1995

1.04

1.12

0.91

0.65

1.23

1.07

1996

0.96

1.04

0.94

0.51

1.20

1.10

1997

1.19

1.19

1.02

0.86

1.38

1.08

1998

1.19

1.27

1.07

0.93

1.45

1.13

1999

1.19

1.27

1.06

0.94

1.46

1.15

2000

0.81

1.04

0.97

0.61

1.22

1.26

2001

0.81

0.88

0.93

0.31

1.07

1.11

2002

1.04

1.19

1.07

0.80

1.39

1.20

2003

1.19

1.19

1.01

0.80

1.39

1.09

2004

1.27

1.35

1.07

0.94

1.53

1.14

2005

1.12

1.27

1.13

0.90

1.48

1.22

2006

1.27

1.35

1.11

0.96

1.50

1.12

2007

1.27

1.35

1.10

0.96

1.47

1.10

2008

1.27

1.27

1.03

0.94

1.47

1.10

2009

1.12

1.35

1.15

0.90

1.47

1.21

2010

1.42

1.50

1.20

0.99

1.69

1.15

2011

1.35

1.50

1.21

1.00

1.67

1.19

2012

1.50

1.50

1.19

0.99

1.69

1.10

2013

1.35

1.42

1.26

0.99

1.71

1.22

2014

1.19

1.35

1.19

0.93

1.58

1.24

2015

1.27

1.35

1.15

0.98

1.53

1.14

2016

1.19

1.27

1.13

0.92

1.51

1.18

2017

1.42

1.50

1.35

1.00

1.76

1.20

2018

1.19

1.35

1.14

0.95

1.51

1.18

2019

1.19

1.42

1.17

0.99

1.59

1.25

Overestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 2, "mm", linf=1.25*65, lmat=26, mk_ratio=2.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

1.12

1.12

0.55

0.01

0.77

0.90

1989

0.81

0.88

0.54

0.01

0.66

0.94

1990

0.81

0.81

0.48

0.00

0.58

0.83

1991

0.96

1.04

0.55

0.01

0.70

0.90

1992

1.04

1.12

0.58

0.01

0.76

0.94

1993

1.04

1.12

0.57

0.01

0.76

0.94

1994

1.12

1.19

0.59

0.02

0.80

0.94

1995

1.04

1.12

0.55

0.00

0.74

0.91

1996

0.96

1.04

0.57

0.01

0.72

0.93

1997

1.19

1.19

0.61

0.03

0.83

0.93

1998

1.19

1.27

0.64

0.06

0.87

0.98

1999

1.19

1.27

0.64

0.06

0.88

0.99

2000

0.81

1.04

0.58

0.02

0.73

1.04

2001

0.81

0.88

0.56

0.01

0.64

0.92

2002

1.04

1.19

0.64

0.06

0.83

1.02

2003

1.19

1.19

0.61

0.02

0.84

0.94

2004

1.27

1.35

0.64

0.07

0.92

0.99

2005

1.12

1.27

0.68

0.11

0.89

1.04

2006

1.27

1.35

0.66

0.08

0.90

0.97

2007

1.27

1.35

0.66

0.06

0.88

0.95

2008

1.27

1.27

0.62

0.04

0.88

0.95

2009

1.12

1.35

0.69

0.11

0.88

1.04

2010

1.42

1.50

0.72

0.21

1.01

1.01

2011

1.35

1.50

0.72

0.25

1.00

1.04

2012

1.50

1.50

0.71

0.20

1.01

0.97

2013

1.35

1.42

0.75

0.32

1.02

1.06

2014

1.19

1.35

0.71

0.22

0.95

1.07

2015

1.27

1.35

0.69

0.09

0.92

0.99

2016

1.19

1.27

0.68

0.16

0.90

1.02

2017

1.42

1.50

0.81

0.32

1.06

1.05

2018

1.19

1.35

0.68

0.12

0.91

1.02

2019

1.19

1.42

0.70

0.19

0.96

1.08

3 Scyliorhinus canicula

Life history parameters:

\(L_{mat}\) = 54.2 cm. \(L_{inf}^{LIT}\) = 75.0 cm. \({M/k}^{LIT}\) = 1.5.

Below, catch/landings at length per year and mean weights-at-length per year are loaded.

wal <- read.csv("Pin8c9awal.csv",sep=",")
head(wal)
##   MeanLength X2014 X2015 X2016 X2017 X2018 X2019
## 1       20.5  23.2  23.2  23.2  23.2  23.2  23.2
## 2       21.5  26.9  26.9  26.9  26.9  26.9  26.9
## 3       22.5  31.0  31.0  31.0  31.0  31.0  31.0
## 4       23.5  35.5  35.5  35.5  35.5  35.5  35.5
## 5       24.5  40.4  40.4  40.4  40.4  40.4  40.4
## 6       25.5  45.7  45.7  45.7  45.7  45.7  45.7
freq <- read.csv("Pin8c9afreq.csv", stringsAsFactors = FALSE)
head(freq)
##   MeanLength X2014 X2015 X2016 X2017 X2018 X2019
## 1       20.5     0     0     0     0     0     0
## 2       21.5     0     0     0     0     0     0
## 3       22.5     0     0     0     0     0     0
## 4       23.5     0     0     0     0     0     0
## 5       24.5     0     0     0     0     0     0
## 6       25.5     0     0     0     0     0     0
bin_plot(freq, 1, "cm")

bin_plot(freq, 2, "cm")

On based on the previous plots of the length distribution BinWidth is fixed equal to 1.

Reference model:

lb_tableSH(freq, 2, "cm", linf=75.0, lmat=54.2, mk_ratio=1.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

2014

0.90

0.94

0.82

0.37

1.07

0.96

2015

0.90

0.90

0.83

0.36

1.08

0.97

2016

0.94

0.94

0.87

0.45

1.10

0.97

2017

0.90

0.90

0.83

0.29

1.07

0.96

2018

0.90

0.90

0.83

0.30

1.07

0.96

2019

0.90

0.90

0.82

0.30

1.06

0.96

Underestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 2, "cm", linf=75.0, lmat=54.2, mk_ratio=0.75*1.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

2014

0.90

0.94

0.82

0.03

0.98

0.94

2015

0.90

0.90

0.83

0.05

0.99

0.94

2016

0.94

0.94

0.87

0.09

1.01

0.94

2017

0.90

0.90

0.83

0.04

0.98

0.94

2018

0.90

0.90

0.83

0.04

0.98

0.94

2019

0.90

0.90

0.82

0.04

0.98

0.93

Overestimated \({M/k}^{LIT}\)

lb_tableSH(freq, 2, "cm", linf=75.0, lmat=54.2, mk_ratio=1.25*1.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

2014

0.90

0.94

0.82

0.78

1.16

0.98

2015

0.90

0.90

0.83

0.72

1.17

0.99

2016

0.94

0.94

0.87

0.84

1.19

0.98

2017

0.90

0.90

0.83

0.66

1.15

0.98

2018

0.90

0.90

0.83

0.67

1.16

0.98

2019

0.90

0.90

0.82

0.70

1.15

0.98

Underestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 2, "cm", linf=0.75*75.0, lmat=54.2, mk_ratio=1.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

2014

0.90

0.94

1.09

0.97

1.43

1.05

2015

0.90

0.90

1.11

0.93

1.44

1.06

2016

0.94

0.94

1.16

0.99

1.47

1.05

2017

0.90

0.90

1.10

0.96

1.42

1.05

2018

0.90

0.90

1.10

0.96

1.42

1.05

2019

0.90

0.90

1.09

0.99

1.42

1.05

Overestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 2, "cm", linf=1.25*75.0, lmat=54.2, mk_ratio=1.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

2014

0.90

0.94

0.65

0.00

0.86

0.89

2015

0.90

0.90

0.66

0.00

0.86

0.89

2016

0.94

0.94

0.70

0.01

0.88

0.89

2017

0.90

0.90

0.66

0.00

0.85

0.89

2018

0.90

0.90

0.66

0.00

0.85

0.89

2019

0.90

0.90

0.65

0.00

0.85

0.88

4 Engraulis encrasicolus

Life history parameters:

\(L_{mat}\) = 11.2 cm \(L_{inf}^{LIT}\) = 18.95 cm. \({M/k}^{LIT}\) = 1.44.

Below, catch/landings at length per year and mean weights-at-length per year are loaded.

MK<-1.44
Linf<-18.95
Lmat<-11.2
wal <- read.csv("wal_anchovy_fv.csv", stringsAsFactors = FALSE)
freq <- read.csv("lengthfreq_anchovy_fv.csv", stringsAsFactors = FALSE)
head(freq)
##   length X1988 X1989    X1990    X1991  X1992   X1993 X1994 X1995     X1996
## 1   3.75     0     0        0        0      0       0     0     0   1349286
## 2   4.25     0     0  4281037   171956   1527   48637     0     0  12676995
## 3   4.75     0     0 18371017  3937228  29018  706893     0     0  67819283
## 4   5.25 64688     0 32251146 54991246  90023 1831793     0     0 160893713
## 5   5.75 86251     0 46583672 80536890 368660 3246934     0     0 129791150
## 6   6.25     0     0 45810220 43303347 983188 5030957     0     0  52811605
##      X1997    X1998    X1999    X2000      X2001    X2002   X2003  X2004  X2005
## 1        0        0        0        0   282386.1    76904       0      0      0
## 2        0        0  1830574   113762   211277.3   274755   35584      0  16257
## 3  1333027  4655928 17055364   806404  1747585.3  1463072  116275  25370 130064
## 4 11492104 25824726 41100084  4886709  5808681.0  3870681  218084  53803 146322
## 5 38721933 57085884 36181025  9251152  9844138.0  8742129  653305 213296  81290
## 6 53185343 82441650 19366091 12753186 12499866.5 13778906 1763017 395900 445378
##    X2006 X2007 X2008 X2009 X2010 X2011  X2012 X2013 X2014 X2015 X2016  X2017
## 1      0     0     0     0     0     0      0     0     0     0     0      0
## 2      0     0     0     0     0     0      0     0     0     0     0      0
## 3      0     0     0     0     0     0      0     0     0     0     0      0
## 4      0     0     0     0     0     0      0     0     0     0     0  45462
## 5      0     0     0     0     0     0      0     0     0     0  9084  34097
## 6 372540     0     0     0     0     0 394206     0     0 80458 54503 217184
##    X2018 X2019
## 1      0     0
## 2      0     0
## 3      0  5047
## 4   5252  5047
## 5 122414     0
## 6 170751 19929
head(wal)
##   length     X1988     X1989     X1990     X1991     X1992     X1993     X1994
## 1   3.75 0.2518403 0.2092267 0.2271301 0.2180922 0.3299293 0.3550483 0.2759662
## 2   4.25 0.3770388 0.3210175 0.3446703 0.3339648 0.4823867 0.5143807 0.4116180
## 3   4.75 0.5396730 0.4696109 0.4993023 0.4877011 0.6760814 0.7150791 0.5872153
## 4   5.25 0.7452009 0.6612985 0.6969456 0.6856975 0.9160444 0.9618145 0.8084307
## 5   5.75 0.9992115 0.9026572 0.9437311 0.9346291 1.2073248 1.2592391 1.0810548
## 6   6.25 1.3074143 1.2005329 1.2459871 1.2414337 1.5549890 1.6119883 1.4109861
##       X1995     X1996     X1997     X1998     X1999     X2000     X2001
## 1 0.2759662 0.2660571 0.2890663 0.2109831 0.2001816 0.2449791 0.2620221
## 2 0.4116180 0.3960943 0.4288719 0.3227909 0.3081657 0.3695502 0.3950267
## 3 0.5872153 0.5641280 0.6089467 0.4710105 0.4521485 0.5325203 0.5689337
## 4 0.8084307 0.7754821 0.8347932 0.6617593 0.6384087 0.7397834 0.7899969
## 5 1.0810548 1.0355821 1.1120032 0.9014165 0.8735282 0.9974136 1.0646578
## 6 1.4109861 1.3499475 1.4462498 1.1966083 1.1643759 1.3116520 1.3995322
##       X2002     X2003     X2004     X2005     X2006     X2007     X2008
## 1 0.2161115 0.2690808 0.2644649 0.2566905 0.2460155 0.2929484 0.3418905
## 2 0.3314068 0.4029565 0.3953278 0.3863321 0.3687180 0.4345656 0.5011364
## 3 0.4845831 0.5769057 0.5650740 0.5555725 0.5282721 0.6169477 0.7039346
## 4 0.6820956 0.7967816 0.7793116 0.7703975 0.7300911 0.8456590 0.9557075
## 5 0.9306895 1.0685792 1.0437756 1.0369636 0.9797234 1.1263528 1.2619092
## 6 1.2373830 1.3984243 1.3643180 1.3615860 1.2828419 1.4647644 1.6280227
##       X2009     X2010     X2011     X2012     X2013     X2014     X2015
## 1 0.2619156 0.3055713 0.2545966 0.2606879 0.2619059 0.2192163 0.2412433
## 2 0.3944381 0.4508095 0.3813615 0.3935505 0.3925561 0.3362669 0.3636509
## 3 0.5675389 0.6368953 0.5461100 0.5674935 0.5624538 0.4918181 0.5236817
## 4 0.7873773 0.8691786 0.7543998 0.7888552 0.7773667 0.6924424 0.7270830
## 5 1.0602915 1.1530702 1.0119242 1.0641714 1.0432063 0.9450092 0.9797736
## 6 1.3927860 1.4940365 1.3245022 1.4001616 1.3660173 1.2566676 1.2878320
##       X2016     X2017     X2018     X2019
## 1 0.2428885 0.3741121 0.3225905 0.3050406
## 2 0.3695925 0.5416183 0.4751011 0.4500820
## 3 0.5367072 0.7524734 0.6701900 0.6359370
## 4 0.7507957 1.0115421 0.9133608 0.8679562
## 5 1.0186707 1.3236670 1.2101712 1.1515516
## 6 1.3473791 1.6936699 1.5662289 1.4921912
bin_plot(freq, 1, "cm")

bin_plot(freq, 2, "cm")

On based on the previous plots of the length distribution BinWidth is fixed equal to 2.

Reference model:

lb_tableSH(freq, 2, "cm", linf=Linf, lmat=Lmat, mk_ratio=MK ,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

0.98

0.98

0.81

0.05

0.93

0.92

1989

0.98

0.98

0.86

0.06

0.93

0.91

1990

0.80

0.62

0.71

0.01

0.84

0.93

1991

0.45

0.62

0.78

0.04

0.76

1.14

1992

0.80

0.80

0.74

0.02

0.88

0.97

1993

1.16

0.80

0.79

0.05

1.03

0.91

1994

0.80

0.80

0.80

0.05

0.88

0.97

1995

0.98

0.98

0.74

0.02

0.91

0.89

1996

0.45

0.45

0.71

0.01

0.53

0.79

1997

0.62

0.62

0.80

0.07

0.79

1.00

1998

0.98

0.80

0.72

0.01

0.91

0.90

1999

0.98

0.62

0.80

0.07

0.98

0.96

2000

0.62

0.80

0.78

0.05

0.80

1.02

2001

0.98

0.98

0.85

0.15

1.00

0.98

2002

0.98

0.80

0.81

0.11

0.97

0.95

2003

0.98

0.98

0.80

0.09

0.95

0.94

2004

0.98

0.98

0.78

0.04

0.94

0.92

2005

0.80

0.80

0.74

0.03

0.86

0.95

2006

0.80

0.80

0.74

0.03

0.87

0.96

2007

0.80

0.80

0.73

0.02

0.86

0.95

2008

0.98

0.98

0.83

0.20

1.01

0.99

2009

0.98

0.98

0.80

0.15

0.98

0.96

2010

0.98

0.98

0.79

0.07

0.94

0.92

2011

0.98

0.98

0.79

0.05

0.93

0.91

2012

0.98

0.98

0.80

0.14

1.00

0.98

2013

0.98

0.98

0.80

0.07

0.95

0.93

2014

0.98

0.98

0.76

0.04

0.92

0.91

2015

0.98

0.98

0.81

0.08

0.93

0.92

2016

0.98

0.98

0.80

0.09

0.95

0.94

2017

0.98

0.98

0.74

0.02

0.89

0.88

2018

0.98

0.98

0.80

0.08

0.96

0.95

2019

0.98

0.98

0.79

0.05

0.95

0.93

Underestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 2, "cm", linf=Linf, lmat=Lmat, mk_ratio=MK*0.75,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

0.98

0.98

0.81

0.01

0.86

0.88

1989

0.98

0.98

0.86

0.03

0.86

0.88

1990

0.80

0.62

0.71

0.00

0.77

0.89

1991

0.45

0.62

0.78

0.01

0.70

1.04

1992

0.80

0.80

0.74

0.00

0.81

0.92

1993

1.16

0.80

0.79

0.00

0.95

0.89

1994

0.80

0.80

0.80

0.00

0.81

0.93

1995

0.98

0.98

0.74

0.00

0.83

0.86

1996

0.45

0.45

0.71

0.00

0.49

0.72

1997

0.62

0.62

0.80

0.00

0.72

0.94

1998

0.98

0.80

0.72

0.00

0.84

0.86

1999

0.98

0.62

0.80

0.00

0.90

0.93

2000

0.62

0.80

0.78

0.00

0.74

0.95

2001

0.98

0.98

0.85

0.03

0.92

0.94

2002

0.98

0.80

0.81

0.01

0.89

0.92

2003

0.98

0.98

0.80

0.01

0.88

0.90

2004

0.98

0.98

0.78

0.00

0.86

0.89

2005

0.80

0.80

0.74

0.00

0.79

0.90

2006

0.80

0.80

0.74

0.00

0.80

0.91

2007

0.80

0.80

0.73

0.00

0.79

0.90

2008

0.98

0.98

0.83

0.02

0.93

0.96

2009

0.98

0.98

0.80

0.00

0.90

0.93

2010

0.98

0.98

0.79

0.00

0.86

0.89

2011

0.98

0.98

0.79

0.00

0.85

0.88

2012

0.98

0.98

0.80

0.01

0.92

0.95

2013

0.98

0.98

0.80

0.00

0.87

0.90

2014

0.98

0.98

0.76

0.00

0.85

0.88

2015

0.98

0.98

0.81

0.01

0.86

0.88

2016

0.98

0.98

0.80

0.01

0.88

0.90

2017

0.98

0.98

0.74

0.00

0.82

0.85

2018

0.98

0.98

0.80

0.00

0.89

0.91

2019

0.98

0.98

0.79

0.00

0.87

0.90

Overestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 2, "cm", linf=Linf, lmat=Lmat, mk_ratio=MK*1.25,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

0.98

0.98

0.81

0.05

1.01

0.94

1989

0.98

0.98

0.86

0.06

1.01

0.94

1990

0.80

0.62

0.71

0.01

0.91

0.97

1991

0.45

0.62

0.78

0.04

0.83

1.22

1992

0.80

0.80

0.74

0.02

0.95

1.01

1993

1.16

0.80

0.79

0.05

1.11

0.92

1994

0.80

0.80

0.80

0.05

0.95

1.01

1995

0.98

0.98

0.74

0.02

0.98

0.91

1996

0.45

0.45

0.71

0.01

0.57

0.84

1997

0.62

0.62

0.80

0.07

0.85

1.05

1998

0.98

0.80

0.72

0.01

0.99

0.92

1999

0.98

0.62

0.80

0.07

1.06

0.99

2000

0.62

0.80

0.78

0.05

0.87

1.07

2001

0.98

0.98

0.85

0.15

1.08

1.00

2002

0.98

0.80

0.81

0.11

1.05

0.98

2003

0.98

0.98

0.80

0.09

1.03

0.96

2004

0.98

0.98

0.78

0.04

1.01

0.94

2005

0.80

0.80

0.74

0.03

0.92

0.98

2006

0.80

0.80

0.74

0.03

0.94

1.00

2007

0.80

0.80

0.73

0.02

0.92

0.98

2008

0.98

0.98

0.83

0.20

1.09

1.02

2009

0.98

0.98

0.80

0.15

1.06

0.98

2010

0.98

0.98

0.79

0.07

1.02

0.95

2011

0.98

0.98

0.79

0.05

1.00

0.93

2012

0.98

0.98

0.80

0.14

1.08

1.00

2013

0.98

0.98

0.80

0.07

1.02

0.95

2014

0.98

0.98

0.76

0.04

1.00

0.93

2015

0.98

0.98

0.81

0.08

1.01

0.94

2016

0.98

0.98

0.80

0.09

1.03

0.96

2017

0.98

0.98

0.74

0.02

0.96

0.90

2018

0.98

0.98

0.80

0.08

1.04

0.97

2019

0.98

0.98

0.79

0.05

1.02

0.95

M/k=1.5.

lb_tableSH(freq, 2, "cm", linf=Linf, lmat=Lmat, mk_ratio=1.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

0.98

0.98

0.81

0.05

0.95

0.92

1989

0.98

0.98

0.86

0.06

0.94

0.92

1990

0.80

0.62

0.71

0.01

0.85

0.94

1991

0.45

0.62

0.78

0.04

0.77

1.15

1992

0.80

0.80

0.74

0.02

0.89

0.98

1993

1.16

0.80

0.79

0.05

1.04

0.91

1994

0.80

0.80

0.80

0.05

0.89

0.98

1995

0.98

0.98

0.74

0.02

0.92

0.89

1996

0.45

0.45

0.71

0.01

0.53

0.80

1997

0.62

0.62

0.80

0.07

0.80

1.01

1998

0.98

0.80

0.72

0.01

0.92

0.90

1999

0.98

0.62

0.80

0.07

1.00

0.97

2000

0.62

0.80

0.78

0.05

0.81

1.03

2001

0.98

0.98

0.85

0.15

1.01

0.98

2002

0.98

0.80

0.81

0.11

0.98

0.96

2003

0.98

0.98

0.80

0.09

0.97

0.94

2004

0.98

0.98

0.78

0.04

0.95

0.92

2005

0.80

0.80

0.74

0.03

0.87

0.95

2006

0.80

0.80

0.74

0.03

0.88

0.97

2007

0.80

0.80

0.73

0.02

0.87

0.95

2008

0.98

0.98

0.83

0.20

1.03

1.00

2009

0.98

0.98

0.80

0.15

0.99

0.97

2010

0.98

0.98

0.79

0.07

0.95

0.93

2011

0.98

0.98

0.79

0.05

0.94

0.91

2012

0.98

0.98

0.80

0.14

1.01

0.98

2013

0.98

0.98

0.80

0.07

0.96

0.93

2014

0.98

0.98

0.76

0.04

0.94

0.91

2015

0.98

0.98

0.81

0.08

0.95

0.92

2016

0.98

0.98

0.80

0.09

0.97

0.94

2017

0.98

0.98

0.74

0.02

0.90

0.88

2018

0.98

0.98

0.80

0.08

0.98

0.95

2019

0.98

0.98

0.79

0.05

0.96

0.93

Underestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 2, "cm", linf=Linf*0.75, lmat=Lmat, mk_ratio=MK,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

0.98

0.98

1.08

0.89

1.24

1.01

1989

0.98

0.98

1.15

0.78

1.24

1.01

1990

0.80

0.62

0.95

0.51

1.12

1.04

1991

0.45

0.62

1.04

0.57

1.02

1.33

1992

0.80

0.80

0.98

0.73

1.17

1.09

1993

1.16

0.80

1.05

0.69

1.37

0.99

1994

0.80

0.80

1.07

0.64

1.17

1.09

1995

0.98

0.98

0.98

0.81

1.21

0.98

1996

0.45

0.45

0.95

0.15

0.70

0.92

1997

0.62

0.62

1.07

0.43

1.05

1.14

1998

0.98

0.80

0.95

0.57

1.22

0.99

1999

0.98

0.62

1.06

0.62

1.31

1.06

2000

0.62

0.80

1.05

0.52

1.07

1.16

2001

0.98

0.98

1.13

0.77

1.33

1.08

2002

0.98

0.80

1.08

0.75

1.29

1.05

2003

0.98

0.98

1.07

0.78

1.27

1.03

2004

0.98

0.98

1.04

0.87

1.25

1.01

2005

0.80

0.80

0.99

0.73

1.14

1.06

2006

0.80

0.80

0.99

0.73

1.16

1.07

2007

0.80

0.80

0.98

0.68

1.14

1.06

2008

0.98

0.98

1.10

0.89

1.35

1.10

2009

0.98

0.98

1.07

0.96

1.31

1.06

2010

0.98

0.98

1.06

0.93

1.25

1.02

2011

0.98

0.98

1.05

0.87

1.24

1.00

2012

0.98

0.98

1.07

0.78

1.33

1.08

2013

0.98

0.98

1.06

0.80

1.26

1.03

2014

0.98

0.98

1.02

0.87

1.23

1.00

2015

0.98

0.98

1.08

0.87

1.24

1.01

2016

0.98

0.98

1.07

0.93

1.27

1.03

2017

0.98

0.98

0.99

0.78

1.19

0.97

2018

0.98

0.98

1.06

0.94

1.29

1.04

2019

0.98

0.98

1.05

0.84

1.26

1.02

Overestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 2, "cm", linf=Linf*1.25, lmat=Lmat, mk_ratio=MK,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1988

0.98

0.98

0.65

0

0.75

0.84

1989

0.98

0.98

0.69

0

0.75

0.84

1990

0.80

0.62

0.57

0

0.67

0.84

1991

0.45

0.62

0.62

0

0.61

1.00

1992

0.80

0.80

0.59

0

0.70

0.88

1993

1.16

0.80

0.63

0

0.82

0.84

1994

0.80

0.80

0.64

0

0.70

0.88

1995

0.98

0.98

0.59

0

0.73

0.81

1996

0.45

0.45

0.57

0

0.42

0.69

1997

0.62

0.62

0.64

0

0.63

0.89

1998

0.98

0.80

0.57

0

0.73

0.82

1999

0.98

0.62

0.64

0

0.79

0.88

2000

0.62

0.80

0.63

0

0.64

0.91

2001

0.98

0.98

0.68

0

0.80

0.89

2002

0.98

0.80

0.65

0

0.78

0.87

2003

0.98

0.98

0.64

0

0.76

0.86

2004

0.98

0.98

0.63

0

0.75

0.84

2005

0.80

0.80

0.59

0

0.68

0.86

2006

0.80

0.80

0.59

0

0.69

0.87

2007

0.80

0.80

0.59

0

0.68

0.86

2008

0.98

0.98

0.66

0

0.81

0.91

2009

0.98

0.98

0.64

0

0.78

0.88

2010

0.98

0.98

0.64

0

0.75

0.84

2011

0.98

0.98

0.63

0

0.74

0.83

2012

0.98

0.98

0.64

0

0.80

0.90

2013

0.98

0.98

0.64

0

0.76

0.85

2014

0.98

0.98

0.61

0

0.74

0.83

2015

0.98

0.98

0.65

0

0.75

0.84

2016

0.98

0.98

0.64

0

0.76

0.86

2017

0.98

0.98

0.59

0

0.71

0.80

2018

0.98

0.98

0.64

0

0.77

0.87

2019

0.98

0.98

0.63

0

0.76

0.85

5 Pollachius pollachius

Life history parameters:

\(L_{mat}\) = 42.3 cm. \(L_{inf}^{LIT}\) = 98.2 cm. \({M/k}^{LIT}\) = 1.648352.

Below, catch/landings at length per year and mean weights-at-length per year are loaded.

freq <- read.csv("Pol89afreq.csv", stringsAsFactors = FALSE)
head(freq)
##   MeanLength X2010 X2011 X2012 X2013 X2014 X2015 X2016  X2017 X2018
## 1       15.5     0     0     0     0     0 0.000     0  1.130     0
## 2       16.5     0     0     0     0     0 0.000     0  0.000     0
## 3       17.5     0     0     0     0     0 0.000     0  2.542     0
## 4       18.5     0     0     0     0     0 0.000     0  3.672     0
## 5       19.5     0     0     0     0     0 4.862     0  6.215     0
## 6       20.5     0     0     0     0     0 4.862     0 12.994     0
wal <- read.csv("Pol89awal.csv", stringsAsFactors = FALSE)
head(wal)
##   MeanLength      X2010      X2011      X2012      X2013      X2014      X2015
## 1       15.5 0.04108271 0.04108271 0.04108271 0.04108271 0.04108271 0.04108271
## 2       16.5 0.04957186 0.04957186 0.04957186 0.04957186 0.04957186 0.04957186
## 3       17.5 0.05915753 0.05915753 0.05915753 0.05915753 0.05915753 0.05915753
## 4       18.5 0.06990645 0.06990645 0.06990645 0.06990645 0.06990645 0.06990645
## 5       19.5 0.08188540 0.08188540 0.08188540 0.08188540 0.08188540 0.08188540
## 6       20.5 0.09516118 0.09516118 0.09516118 0.09516118 0.09516118 0.09516118
##        X2016      X2017      X2018
## 1 0.04108271 0.04108271 0.04108271
## 2 0.04957186 0.04957186 0.04957186
## 3 0.05915753 0.05915753 0.05915753
## 4 0.06990645 0.06990645 0.06990645
## 5 0.08188540 0.08188540 0.08188540
## 6 0.09516118 0.09516118 0.09516118
bin_plot(freq, 4, "cm")

On based on the previous plots of the length distribution BinWidth is fixed equal to 4.

Reference model:

lb_tableSH(freq, 4, "cm", linf=98.2, lmat=42.3, mk_ratio=1.648352,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

2010

0.80

0.80

0.72

0.03

0.67

0.87

2011

1.18

1.18

0.80

0.07

0.91

0.94

2012

1.18

1.18

0.80

0.11

0.95

0.98

2013

1.09

1.09

0.79

0.12

0.89

0.97

2014

0.80

0.99

0.77

0.10

0.83

1.08

2015

0.71

0.80

0.74

0.04

0.70

0.96

2016

0.80

0.90

0.71

0.03

0.71

0.93

2017

0.80

0.90

0.80

0.09

0.79

1.02

2018

0.80

0.99

0.82

0.11

0.86

1.12

Underestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 4, "cm", linf=98.2, lmat=42.3, mk_ratio=0.75*1.648352,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

2010

0.80

0.80

0.72

0.01

0.61

0.81

2011

1.18

1.18

0.80

0.02

0.83

0.90

2012

1.18

1.18

0.80

0.03

0.86

0.94

2013

1.09

1.09

0.79

0.01

0.81

0.92

2014

0.80

0.99

0.77

0.01

0.76

1.01

2015

0.71

0.80

0.74

0.01

0.64

0.89

2016

0.80

0.90

0.71

0.01

0.65

0.86

2017

0.80

0.90

0.80

0.03

0.72

0.95

2018

0.80

0.99

0.82

0.05

0.79

1.04

Overestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 4, "cm", linf=98.2, lmat=42.3, mk_ratio=1.25*1.648352,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

2010

0.80

0.80

0.72

0.07

0.73

0.91

2011

1.18

1.18

0.80

0.11

0.99

0.97

2012

1.18

1.18

0.80

0.22

1.03

1.01

2013

1.09

1.09

0.79

0.18

0.97

1.00

2014

0.80

0.99

0.77

0.16

0.91

1.14

2015

0.71

0.80

0.74

0.07

0.76

1.02

2016

0.80

0.90

0.71

0.04

0.78

0.97

2017

0.80

0.90

0.80

0.14

0.86

1.08

2018

0.80

0.99

0.82

0.19

0.94

1.17

M/k=1.5.

lb_tableSH(freq, 4, "cm", linf=98.2, lmat=42.3, mk_ratio=1.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

2010

0.80

0.80

0.72

0.02

0.65

0.85

2011

1.18

1.18

0.80

0.04

0.88

0.93

2012

1.18

1.18

0.80

0.06

0.91

0.97

2013

1.09

1.09

0.79

0.08

0.86

0.95

2014

0.80

0.99

0.77

0.05

0.81

1.06

2015

0.71

0.80

0.74

0.02

0.68

0.94

2016

0.80

0.90

0.71

0.02

0.69

0.90

2017

0.80

0.90

0.80

0.06

0.76

1.00

2018

0.80

0.99

0.82

0.07

0.83

1.09

Underestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 4, "cm", linf=0.75*98.2, lmat=42.3, mk_ratio=1.648352,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

2010

0.80

0.80

0.96

0.14

0.89

0.98

2011

1.18

1.18

1.06

0.66

1.21

1.04

2012

1.18

1.18

1.07

0.69

1.26

1.08

2013

1.09

1.09

1.05

0.52

1.19

1.07

2014

0.80

0.99

1.02

0.58

1.11

1.22

2015

0.71

0.80

0.98

0.23

0.93

1.10

2016

0.80

0.90

0.94

0.18

0.95

1.05

2017

0.80

0.90

1.07

0.33

1.05

1.16

2018

0.80

0.99

1.09

0.61

1.15

1.26

Overestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 4, "cm", linf=1.25*98.2, lmat=42.3, mk_ratio=1.648352,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

2010

0.80

0.80

0.58

0.00

0.54

0.78

2011

1.18

1.18

0.64

0.01

0.73

0.86

2012

1.18

1.18

0.64

0.00

0.76

0.89

2013

1.09

1.09

0.63

0.00

0.71

0.88

2014

0.80

0.99

0.61

0.00

0.67

0.97

2015

0.71

0.80

0.59

0.00

0.56

0.86

2016

0.80

0.90

0.57

0.00

0.57

0.83

2017

0.80

0.90

0.64

0.00

0.63

0.92

2018

0.80

0.99

0.65

0.00

0.69

1.00

6 Pagellus bogaraveo

Life history parameters:

\(L_{mat}\) = 33 cm. \(L_{inf}^{LIT}\) = 62 cm. \({M/k}^{LIT}\) = 1.42.

Below, catch/landings at length per year and mean weights-at-length per year are loaded.

wal <- read.csv("SBR9_mw.csv", stringsAsFactors = FALSE)
freq <- read.csv("SBR9_numbers.csv", stringsAsFactors = FALSE)
Linf<-62
Lmat<-33
MK<-1.42

head(freq)
##   MeanLength X1997 X1998 X1999 X2000 X2001 X2002 X2003 X2004 X2005 X2006 X2007
## 1       24.5   422   170     0     0     0     0     0     0     0     0     0
## 2       25.5  1553   718     0     0     0     0     0     0     0     0     0
## 3       26.5  9222  6454   527    48    42    91    52     0     0     0     0
## 4       27.5 33437 25121  2840   717   912   376   599     0  1012  1583   138
## 5       28.5 69291 53905  9521  3875  3214  2266  3516   254  7417 11301  1224
## 6       29.5 81193 66854 13037 11083 11024  8050 10162  2416 15187 27559  9655
##   X2008 X2009 X2010 X2011 X2012 X2013 X2014 X2015 X2016 X2017 X2018 X2019
## 1     0     0     0     0     0     0     0     0     0     0     0     0
## 2     0     0     0     0     0     0     0     0     0     0     0     0
## 3     0     0     0     0     0     0     0     0     0     0     0     0
## 4     0     0     0     0     0     0     0     0     0     0     0     0
## 5   160   267   285   742     0     0     0     0     0     0     0     0
## 6  6810  2139  7243 12657   872     0     0     0   160     0     0     0
head(wal)
##   MeanLength X1997 X1998 X1999 X2000 X2001 X2002 X2003 X2004 X2005 X2006 X2007
## 1       24.5   200   200   200   200   200   200   200   200   200   200   200
## 2       25.5   227   227   227   227   227   227   227   227   227   227   227
## 3       26.5   256   256   256   256   256   256   256   256   256   256   256
## 4       27.5   288   288   288   288   288   288   288   288   288   288   288
## 5       28.5   322   322   322   322   322   322   322   322   322   322   322
## 6       29.5   359   359   359   359   359   359   359   359   359   359   359
##   X2008 X2009 X2010 X2011 X2012 X2013 X2014 X2015 X2016 X2017 X2018 X2019
## 1   200   200   200   200   200   200   200   200   200   200   200   200
## 2   227   227   227   227   227   227   227   227   227   227   227   227
## 3   256   256   256   256   256   256   256   256   256   256   256   256
## 4   288   288   288   288   288   288   288   288   288   288   288   288
## 5   322   322   322   322   322   322   322   322   322   322   322   322
## 6   359   359   359   359   359   359   359   359   359   359   359   359
bin_plot(freq, 1, "cm")

bin_plot(freq, 2, "cm")

On based on the previous plots of the length distribution BinWidth is fixed equal to 2.

Reference model:

lb_tableSH(freq, 2, "cm", linf=Linf, lmat=Lmat, mk_ratio=MK ,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1997

0.88

0.94

0.82

0.11

0.86

0.97

1998

0.88

0.94

0.78

0.04

0.82

0.92

1999

0.94

1.00

0.80

0.07

0.87

0.94

2000

1.00

1.00

0.80

0.06

0.90

0.93

2001

0.94

1.00

0.82

0.09

0.89

0.96

2002

0.94

1.00

0.83

0.11

0.92

0.99

2003

0.94

1.00

0.83

0.13

0.93

1.00

2004

0.94

1.00

0.82

0.08

0.87

0.94

2005

0.94

1.00

0.82

0.11

0.89

0.95

2006

0.94

0.94

0.81

0.07

0.87

0.93

2007

0.94

1.00

0.83

0.11

0.89

0.96

2008

0.94

1.00

0.84

0.14

0.90

0.97

2009

0.94

1.00

0.85

0.16

0.91

0.98

2010

0.94

1.00

0.82

0.08

0.86

0.93

2011

0.94

0.94

0.84

0.12

0.87

0.94

2012

0.94

1.00

0.84

0.11

0.87

0.93

2013

1.00

1.00

0.75

0.02

0.84

0.87

2014

1.00

1.00

0.83

0.08

0.89

0.93

2015

1.00

1.06

0.85

0.15

0.94

0.97

2016

0.94

1.00

0.82

0.09

0.89

0.96

2017

1.00

1.06

0.80

0.06

0.90

0.93

2018

1.00

1.06

0.81

0.07

0.90

0.93

2019

1.00

1.00

0.79

0.05

0.89

0.92

Underestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 2, "cm", linf=Linf, lmat=Lmat, mk_ratio=MK*0.75,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1997

0.88

0.94

0.82

0.03

0.79

0.92

1998

0.88

0.94

0.78

0.01

0.76

0.88

1999

0.94

1.00

0.80

0.02

0.80

0.90

2000

1.00

1.00

0.80

0.02

0.82

0.89

2001

0.94

1.00

0.82

0.03

0.82

0.92

2002

0.94

1.00

0.83

0.04

0.85

0.95

2003

0.94

1.00

0.83

0.04

0.85

0.95

2004

0.94

1.00

0.82

0.03

0.80

0.90

2005

0.94

1.00

0.82

0.03

0.82

0.91

2006

0.94

0.94

0.81

0.02

0.80

0.89

2007

0.94

1.00

0.83

0.04

0.82

0.92

2008

0.94

1.00

0.84

0.05

0.83

0.93

2009

0.94

1.00

0.85

0.06

0.84

0.94

2010

0.94

1.00

0.82

0.03

0.79

0.88

2011

0.94

0.94

0.84

0.05

0.80

0.90

2012

0.94

1.00

0.84

0.05

0.80

0.89

2013

1.00

1.00

0.75

0.01

0.77

0.84

2014

1.00

1.00

0.83

0.04

0.82

0.89

2015

1.00

1.06

0.85

0.07

0.86

0.93

2016

0.94

1.00

0.82

0.03

0.82

0.92

2017

1.00

1.06

0.80

0.02

0.83

0.89

2018

1.00

1.06

0.81

0.02

0.83

0.90

2019

1.00

1.00

0.79

0.02

0.81

0.88

Overestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 2, "cm", linf=Linf, lmat=Lmat, mk_ratio=MK*1.25,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1997

0.88

0.94

0.82

0.19

0.93

1.00

1998

0.88

0.94

0.78

0.10

0.89

0.96

1999

0.94

1.00

0.80

0.14

0.94

0.97

2000

1.00

1.00

0.80

0.14

0.97

0.96

2001

0.94

1.00

0.82

0.18

0.97

0.99

2002

0.94

1.00

0.83

0.24

0.99

1.02

2003

0.94

1.00

0.83

0.30

1.00

1.03

2004

0.94

1.00

0.82

0.20

0.94

0.97

2005

0.94

1.00

0.82

0.22

0.96

0.99

2006

0.94

0.94

0.81

0.16

0.94

0.97

2007

0.94

1.00

0.83

0.20

0.96

0.99

2008

0.94

1.00

0.84

0.24

0.97

1.00

2009

0.94

1.00

0.85

0.25

0.98

1.01

2010

0.94

1.00

0.82

0.13

0.93

0.96

2011

0.94

0.94

0.84

0.18

0.94

0.97

2012

0.94

1.00

0.84

0.17

0.94

0.96

2013

1.00

1.00

0.75

0.05

0.91

0.90

2014

1.00

1.00

0.83

0.15

0.96

0.95

2015

1.00

1.06

0.85

0.28

1.01

1.00

2016

0.94

1.00

0.82

0.19

0.96

0.99

2017

1.00

1.06

0.80

0.16

0.97

0.96

2018

1.00

1.06

0.81

0.18

0.97

0.96

2019

1.00

1.00

0.79

0.13

0.96

0.95

M/k=1.5.

lb_tableSH(freq, 2, "cm", linf=Linf, lmat=Lmat, mk_ratio=1.5,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1997

0.88

0.94

0.82

0.11

0.88

0.98

1998

0.88

0.94

0.78

0.04

0.84

0.93

1999

0.94

1.00

0.80

0.07

0.89

0.95

2000

1.00

1.00

0.80

0.06

0.91

0.94

2001

0.94

1.00

0.82

0.09

0.91

0.97

2002

0.94

1.00

0.83

0.11

0.94

1.00

2003

0.94

1.00

0.83

0.13

0.94

1.01

2004

0.94

1.00

0.82

0.08

0.89

0.95

2005

0.94

1.00

0.82

0.11

0.90

0.96

2006

0.94

0.94

0.81

0.07

0.88

0.94

2007

0.94

1.00

0.83

0.11

0.91

0.97

2008

0.94

1.00

0.84

0.14

0.92

0.98

2009

0.94

1.00

0.85

0.16

0.93

0.99

2010

0.94

1.00

0.82

0.08

0.87

0.93

2011

0.94

0.94

0.84

0.12

0.89

0.95

2012

0.94

1.00

0.84

0.11

0.88

0.94

2013

1.00

1.00

0.75

0.02

0.86

0.88

2014

1.00

1.00

0.83

0.08

0.91

0.93

2015

1.00

1.06

0.85

0.15

0.95

0.98

2016

0.94

1.00

0.82

0.09

0.91

0.97

2017

1.00

1.06

0.80

0.06

0.91

0.94

2018

1.00

1.06

0.81

0.07

0.92

0.94

2019

1.00

1.00

0.79

0.05

0.90

0.93

Underestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 2, "cm", linf=Linf*0.75, lmat=Lmat, mk_ratio=MK,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1997

0.88

0.94

1.09

0.52

1.15

1.08

1998

0.88

0.94

1.04

0.42

1.10

1.03

1999

0.94

1.00

1.07

0.61

1.17

1.05

2000

1.00

1.00

1.07

0.72

1.19

1.03

2001

0.94

1.00

1.09

0.70

1.19

1.07

2002

0.94

1.00

1.11

0.72

1.23

1.10

2003

0.94

1.00

1.11

0.69

1.24

1.11

2004

0.94

1.00

1.09

0.52

1.16

1.05

2005

0.94

1.00

1.10

0.57

1.18

1.06

2006

0.94

0.94

1.08

0.54

1.16

1.04

2007

0.94

1.00

1.11

0.62

1.19

1.07

2008

0.94

1.00

1.12

0.64

1.20

1.08

2009

0.94

1.00

1.13

0.69

1.21

1.09

2010

0.94

1.00

1.10

0.57

1.15

1.03

2011

0.94

0.94

1.12

0.50

1.16

1.05

2012

0.94

1.00

1.12

0.51

1.15

1.04

2013

1.00

1.00

1.00

0.45

1.12

0.97

2014

1.00

1.00

1.11

0.67

1.19

1.03

2015

1.00

1.06

1.13

0.79

1.25

1.08

2016

0.94

1.00

1.09

0.69

1.19

1.07

2017

1.00

1.06

1.07

0.80

1.20

1.04

2018

1.00

1.06

1.08

0.78

1.20

1.04

2019

1.00

1.00

1.06

0.74

1.18

1.02

Overestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 2, "cm", linf=Linf*1.25, lmat=Lmat, mk_ratio=MK,wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1997

0.88

0.94

0.65

0

0.69

0.87

1998

0.88

0.94

0.62

0

0.66

0.83

1999

0.94

1.00

0.64

0

0.70

0.85

2000

1.00

1.00

0.64

0

0.72

0.85

2001

0.94

1.00

0.65

0

0.71

0.87

2002

0.94

1.00

0.67

0

0.74

0.90

2003

0.94

1.00

0.67

0

0.74

0.90

2004

0.94

1.00

0.65

0

0.70

0.85

2005

0.94

1.00

0.66

0

0.71

0.87

2006

0.94

0.94

0.65

0

0.69

0.85

2007

0.94

1.00

0.67

0

0.71

0.87

2008

0.94

1.00

0.67

0

0.72

0.88

2009

0.94

1.00

0.68

0

0.73

0.89

2010

0.94

1.00

0.66

0

0.69

0.84

2011

0.94

0.94

0.67

0

0.70

0.85

2012

0.94

1.00

0.67

0

0.69

0.85

2013

1.00

1.00

0.60

0

0.67

0.80

2014

1.00

1.00

0.67

0

0.71

0.84

2015

1.00

1.06

0.68

0

0.75

0.88

2016

0.94

1.00

0.65

0

0.71

0.87

2017

1.00

1.06

0.64

0

0.72

0.85

2018

1.00

1.06

0.65

0

0.72

0.85

2019

1.00

1.00

0.64

0

0.71

0.84

7 Trisopterus luscus

Life history parameters:

\(L_{mat}\) = 19.2 cm. \(L_{inf}^{LIT}\) = 46.7 cm. \({M/k}^{LIT}\) = 1.761846.

Below, catch/landings at length per year and mean weights-at-length per year are loaded.

LHTs <- data.frame(t(c("Trisopterus luscus", 19.2, 1, 46.7, 0.21, 1.27, 9)))
colnames(LHTs) <- c("Sp", "L50", "A50", "Linf", "K", "t0", "Amax")

LHTs$MKvbgf <- 1.761846
LHTs$MKnls <- 3.111105
LHTs$MK <- 1.5

library(readr)
faneca <- read_csv("faneca_ieo.csv")
sizeDat <- na.omit(faneca %>%
                    select(TLength, Year))
anho <- as.character(unique(sort(faneca$Year)))
sizeDat <- data.frame(sizeDat)
br = seq(min(sizeDat$TLength), max(sizeDat$TLength) + 1, by = 1)
ranges = paste(head(br, -1), br[-1], sep = " - ")

sizeDat <- unstack(sizeDat)

freq <- data.frame(matrix(NA, nrow = length(br)-1, ncol = length(anho)+1))

for (i in 1:length(anho)) {
  freq[, 1] <- head(br, -1)
  freq[, i + 1] <-
    hist(
      sizeDat[[i]],
      breaks = br - 1,
      include.lowest = TRUE,
      plot = FALSE
    )$counts
  colnames(freq) <- c("MeanLength", anho)
}

freq[,1]=seq(10.5,43.5)
br=seq(10.5,43.5)
head(freq)
##   MeanLength 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011
## 1       10.5    0    0    0    0    0    1    0    0    0    0    0    0    0
## 2       11.5    0    0    0    0    0    1    0    0    0    0    0    0    0
## 3       12.5    0    0    1    0    0    0    0    0    0    0    0    0    0
## 4       13.5    0    0    0    2    1    1    0    2    4    2    0    0    1
## 5       14.5    0    2    5    5    2    2    1    2   12    1    1    1    0
## 6       15.5    0    0    8    3   14    2   10    2   18    0    3    1    6
##   2012 2013 2014 2015 2016 2017 2018
## 1    0    0    1    0    0    0    0
## 2    0    0    0    0    0    2    0
## 3    0    0    0    0    0    1    0
## 4    0    0    1    0    0    0    0
## 5    0    2    0    0    0    3    2
## 6    1    3    5    9    3    3    5
alfaWL=0.01303439
betaWL=2.969518

wal <- freq
for (i in 2:dim(freq)[2]) {
  wal[, i] <- alfaWL * ((br) ^ betaWL)
}
colnames(wal) <- c("MeanLength", anho)
head(wal)
##   MeanLength     1999     2000     2001     2002     2003     2004     2005
## 1       10.5 14.04529 14.04529 14.04529 14.04529 14.04529 14.04529 14.04529
## 2       11.5 18.40145 18.40145 18.40145 18.40145 18.40145 18.40145 18.40145
## 3       12.5 23.57137 23.57137 23.57137 23.57137 23.57137 23.57137 23.57137
## 4       13.5 29.62356 29.62356 29.62356 29.62356 29.62356 29.62356 29.62356
## 5       14.5 36.62637 36.62637 36.62637 36.62637 36.62637 36.62637 36.62637
## 6       15.5 44.64801 44.64801 44.64801 44.64801 44.64801 44.64801 44.64801
##       2006     2007     2008     2009     2010     2011     2012     2013
## 1 14.04529 14.04529 14.04529 14.04529 14.04529 14.04529 14.04529 14.04529
## 2 18.40145 18.40145 18.40145 18.40145 18.40145 18.40145 18.40145 18.40145
## 3 23.57137 23.57137 23.57137 23.57137 23.57137 23.57137 23.57137 23.57137
## 4 29.62356 29.62356 29.62356 29.62356 29.62356 29.62356 29.62356 29.62356
## 5 36.62637 36.62637 36.62637 36.62637 36.62637 36.62637 36.62637 36.62637
## 6 44.64801 44.64801 44.64801 44.64801 44.64801 44.64801 44.64801 44.64801
##       2014     2015     2016     2017     2018
## 1 14.04529 14.04529 14.04529 14.04529 14.04529
## 2 18.40145 18.40145 18.40145 18.40145 18.40145
## 3 23.57137 23.57137 23.57137 23.57137 23.57137
## 4 29.62356 29.62356 29.62356 29.62356 29.62356
## 5 36.62637 36.62637 36.62637 36.62637 36.62637
## 6 44.64801 44.64801 44.64801 44.64801 44.64801
bin_plot(freq, 1, "cm")

On based on the previous plots of the length distribution BinWidth is fixed equal to 1.

Reference model:

lb_tableSH(freq, 1, "cm", 46.7 , 19.2, 1.761846, wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1999

1.07

1.12

0.65

0.00

0.79

0.89

2000

0.96

1.02

0.60

0.00

0.74

0.88

2001

0.91

1.02

0.64

0.01

0.75

0.92

2002

0.96

1.02

0.65

0.01

0.77

0.92

2003

0.96

1.07

0.65

0.01

0.78

0.92

2004

1.02

1.07

0.62

0.00

0.78

0.90

2005

0.96

1.02

0.63

0.00

0.76

0.90

2006

0.96

1.02

0.63

0.00

0.75

0.89

2007

0.96

1.02

0.62

0.00

0.75

0.89

2008

1.02

1.02

0.60

0.00

0.75

0.87

2009

1.02

1.07

0.63

0.00

0.78

0.89

2010

1.02

1.02

0.66

0.01

0.75

0.87

2011

0.96

1.02

0.63

0.01

0.75

0.89

2012

1.02

1.07

0.63

0.00

0.76

0.88

2013

0.96

1.02

0.58

0.00

0.73

0.87

2014

1.02

1.07

0.59

0.00

0.75

0.87

2015

1.02

1.07

0.68

0.02

0.80

0.92

2016

1.02

1.07

0.63

0.01

0.76

0.88

2017

1.02

1.07

0.63

0.00

0.78

0.90

2018

1.02

1.07

0.63

0.00

0.78

0.89

Underestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 1,"cm",46.7 , 19.2, 0.75*1.761846, wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1999

1.07

1.12

0.65

0

0.72

0.84

2000

0.96

1.02

0.60

0

0.67

0.83

2001

0.91

1.02

0.64

0

0.68

0.86

2002

0.96

1.02

0.65

0

0.70

0.87

2003

0.96

1.07

0.65

0

0.70

0.87

2004

1.02

1.07

0.62

0

0.70

0.85

2005

0.96

1.02

0.63

0

0.69

0.85

2006

0.96

1.02

0.63

0

0.68

0.84

2007

0.96

1.02

0.62

0

0.68

0.84

2008

1.02

1.02

0.60

0

0.68

0.82

2009

1.02

1.07

0.63

0

0.70

0.85

2010

1.02

1.02

0.66

0

0.68

0.82

2011

0.96

1.02

0.63

0

0.68

0.84

2012

1.02

1.07

0.63

0

0.69

0.83

2013

0.96

1.02

0.58

0

0.66

0.82

2014

1.02

1.07

0.59

0

0.68

0.82

2015

1.02

1.07

0.68

0

0.73

0.87

2016

1.02

1.07

0.63

0

0.69

0.83

2017

1.02

1.07

0.63

0

0.71

0.85

2018

1.02

1.07

0.63

0

0.70

0.85

Overestimated \({M/k}^{LIT}\).

lb_tableSH(freq, 1,"cm",46.7 , 19.2, 1.25*1.761846, wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1999

1.07

1.12

0.65

0.03

0.87

0.92

2000

0.96

1.02

0.60

0.01

0.81

0.92

2001

0.91

1.02

0.64

0.02

0.82

0.96

2002

0.96

1.02

0.65

0.02

0.84

0.96

2003

0.96

1.07

0.65

0.02

0.85

0.96

2004

1.02

1.07

0.62

0.01

0.85

0.93

2005

0.96

1.02

0.63

0.01

0.83

0.94

2006

0.96

1.02

0.63

0.02

0.82

0.93

2007

0.96

1.02

0.62

0.01

0.82

0.93

2008

1.02

1.02

0.60

0.01

0.82

0.90

2009

1.02

1.07

0.63

0.01

0.85

0.93

2010

1.02

1.02

0.66

0.03

0.82

0.90

2011

0.96

1.02

0.63

0.01

0.82

0.93

2012

1.02

1.07

0.63

0.01

0.84

0.92

2013

0.96

1.02

0.58

0.00

0.80

0.90

2014

1.02

1.07

0.59

0.01

0.82

0.90

2015

1.02

1.07

0.68

0.04

0.87

0.96

2016

1.02

1.07

0.63

0.01

0.83

0.91

2017

1.02

1.07

0.63

0.02

0.85

0.94

2018

1.02

1.07

0.63

0.02

0.85

0.93

M/k=1.5.

lb_tableSH(freq, 1,"cm",46.7 , 19.2, 1.5, wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1999

1.07

1.12

0.65

0

0.75

0.86

2000

0.96

1.02

0.60

0

0.70

0.85

2001

0.91

1.02

0.64

0

0.71

0.89

2002

0.96

1.02

0.65

0

0.73

0.89

2003

0.96

1.07

0.65

0

0.73

0.89

2004

1.02

1.07

0.62

0

0.73

0.87

2005

0.96

1.02

0.63

0

0.72

0.87

2006

0.96

1.02

0.63

0

0.71

0.86

2007

0.96

1.02

0.62

0

0.71

0.86

2008

1.02

1.02

0.60

0

0.71

0.84

2009

1.02

1.07

0.63

0

0.73

0.87

2010

1.02

1.02

0.66

0

0.71

0.84

2011

0.96

1.02

0.63

0

0.71

0.86

2012

1.02

1.07

0.63

0

0.72

0.86

2013

0.96

1.02

0.58

0

0.69

0.84

2014

1.02

1.07

0.59

0

0.71

0.84

2015

1.02

1.07

0.68

0

0.76

0.90

2016

1.02

1.07

0.63

0

0.72

0.85

2017

1.02

1.07

0.63

0

0.74

0.87

2018

1.02

1.07

0.63

0

0.73

0.87

Underestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 1,"cm",0.75*46.7 , 19.2, 1.761846, wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1999

1.07

1.12

0.87

0.27

1.06

0.99

2000

0.96

1.02

0.80

0.14

0.99

0.98

2001

0.91

1.02

0.86

0.23

1.00

1.03

2002

0.96

1.02

0.86

0.27

1.03

1.03

2003

0.96

1.07

0.87

0.29

1.04

1.03

2004

1.02

1.07

0.82

0.28

1.04

1.00

2005

0.96

1.02

0.84

0.22

1.01

1.01

2006

0.96

1.02

0.84

0.19

1.00

0.99

2007

0.96

1.02

0.82

0.18

1.00

0.99

2008

1.02

1.02

0.80

0.16

1.00

0.96

2009

1.02

1.07

0.84

0.25

1.03

0.99

2010

1.02

1.02

0.88

0.15

1.00

0.97

2011

0.96

1.02

0.84

0.15

1.00

0.99

2012

1.02

1.07

0.84

0.21

1.02

0.98

2013

0.96

1.02

0.77

0.13

0.97

0.97

2014

1.02

1.07

0.79

0.15

1.00

0.97

2015

1.02

1.07

0.90

0.37

1.07

1.03

2016

1.02

1.07

0.84

0.20

1.02

0.98

2017

1.02

1.07

0.85

0.27

1.04

1.00

2018

1.02

1.07

0.84

0.24

1.03

1.00

Overestimated \(L_{inf}^{LIT}\).

lb_tableSH(freq, 1,"cm",1.25*46.7 , 19.2, 1.761846, wal)

Conservation

Optimizing Yield

MSY

Year

Lc / Lmat

L25% / Lmat

Lmax 5 / Linf

Pmega

Lmean / Lopt

Lmean / LF = M

1999

1.07

1.12

0.52

0

0.64

0.81

2000

0.96

1.02

0.48

0

0.59

0.80

2001

0.91

1.02

0.52

0

0.60

0.83

2002

0.96

1.02

0.52

0

0.62

0.83

2003

0.96

1.07

0.52

0

0.62

0.84

2004

1.02

1.07

0.49

0

0.62

0.81

2005

0.96

1.02

0.50

0

0.61

0.82

2006

0.96

1.02

0.50

0

0.60

0.81

2007

0.96

1.02

0.49

0

0.60

0.81

2008

1.02

1.02

0.48

0

0.60

0.79

2009

1.02

1.07

0.50

0

0.62

0.81

2010

1.02

1.02

0.53

0

0.60

0.79

2011

0.96

1.02

0.51

0

0.60

0.80

2012

1.02

1.07

0.50

0

0.61

0.80

2013

0.96

1.02

0.46

0

0.58

0.79

2014

1.02

1.07

0.47

0

0.60

0.79

2015

1.02

1.07

0.54

0

0.64

0.84

2016

1.02

1.07

0.51

0

0.61

0.80

2017

1.02

1.07

0.51

0

0.63

0.82

2018

1.02

1.07

0.51

0

0.62

0.81